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Abstract 

Two-mode  nonlinear  interaction  (second-harmonic  and 
second-subharmonic  generation)  in  a  planar  waveguide 
with  a  small  periodic  corrugation  at  the  surface  is  stud¬ 
ied.  Scattering  of  the  interacting  fields  on  the  corruga¬ 
tion  leads  to  field  localization  that  enhances  the  nonlin¬ 
ear  process  provided  that  all  the  interactions  are  phase 
matched.  Conditions  for  the  overall  phase  matching  are 
found.  Compared  with  a  perfectly  quasi-phase  matched 
waveguide,  better  values  of  squeezing  as  well  as  higher 
intensities  are  reached  under  these  conditions.  In  pulsed 
regime,  eigenmodes  suitable  for  squeezed-light  generation 
have  been  found. 


I.  INTRODUCTION 

Since  the  pioneering  work  by  Armstrong  [1]  on  the  pro¬ 
cess  of  second-harmonic  generation  has  occurred,  spatio- 
temporal  properties  of  the  nonlinearly  interacting  clas¬ 
sical  fields  have  been  studied  in  detail  by  many  authors 
both  theoretically  and  experimentally.  A  new  impulse  in 
these  studies  has  occurred  when  people  understood  that 
this  process  can  give  rise  to  the  fields  with  nonclassical 
properties  (for  a  review,  see,  e.g.  [2-5]).  Namely  light 
with  electric-field  amplitude  fluctuations  suppressed  be¬ 
low  the  limit  given  by  quantum  mechanics  [6]  can  be  gen¬ 
erated  both  in  the  pump  and  second-subharmonic  fields. 
Also  light  with  nonclassical  photon-number  statistics  can 
be  obtained  -  pairwise  character  of  photon-number  statis¬ 
tics  [7,  8]  generated  in  the  spontaneous  process  of  second- 
subharmonic  generation  has  been  observed  [9]. 

It  has  been  shown  that  the  best  conditions  for 
squeezed-light  generation  in  homogeneous  nonlinear  me¬ 
dia  occur,  provided  that  the  nonlinear  two-mode  interac¬ 
tion  is  perfectly  phase  matched.  Under  these  conditions, 
the  principal  squeeze  variance  of  the  second-subharmonic 
field  can  asymptotically  reach  zero  when  the  gain  of  the 
nonlinear  interaction  increases.  On  the  other  hand,  the 
pump-field  principal  squeeze  variance  cannot  be  less  than 
0.5  [10].  If  large  values  of  the  nonlinear  phase  mismatch 
are  allowed,  this  limit  can  be  overcome  due  to  a  nonlinear 
phase  modulation,  as  suggested  in  [11,  12].  However,  the 
generated  signal  is  very  weak. 

The  most  common  method  how  to  compensate  for 
the  natural  nonlinear  phase  mismatch  that  occurs  in 
commonly  used  nonlinear  materials  is  to  introduce  an 
additional  periodic  modulation  of  the  susceptibil¬ 
ity  using  periodical  poling  [13-15].  Several  methods 
for  the  additional  modulation  of  the  local  amplitude  of 
this  quasi-phase-matched  interaction  have  been  devel¬ 
oped  [16].  These  methods  allow  to  reach  a  spectrally 
broad-band  two-mode  interaction  and  so  femtosecond 
pumping  of  the  nonlinear  process  is  possible  [17]. 

In  order  to  effectively  increase  low  values  of  the  nonlin¬ 
ear  interaction  in  real  materials,  configurations  in  which 
a  nonlinear  medium  is  put  inside  a  cavity  are  usually  used 
to  generate  squeezed  light  (e.g.  [6,  18]). 


In  a  waveguiding  geometry  that  profits  from  a  strong 
spatial  localization  of  the  interacting  optical  fields  re¬ 
sulting  in  high  values  of  the  effective  nonlinearity,  an¬ 
other  method  to  reach  a  nonlinear  phase  mismatch  is 
possible.  One  of  the  nonlinearly  interacting  fields  can  be 
coupled  through  its  evanescent  waves  into  another  field  of 
the  same  frequency  propagating  in  a  neighbouring  waveg¬ 
uide.  An  exchange  of  energy  between  these  two  linearly 
coupled  fields  introduces  a  spatial  modulation  of  the  non¬ 
linearly  interacting  field  that  can  be  set  to  compensate  for 
the  nonlinear  phase  mismatch  [19].  Interaction  of  fields  in 
different  waveguides  through  their  evanescent  waves  can 
be  used  in  various  configurations  that  modify  nonclassi¬ 
cal  properties  of  optical  fields  emerging  from  nonlinear 
interactions  [20]. 

A  waveguiding  geometry  allows  another  possibility  to 
tailor  the  nonlinear  process  -  a  linear  periodic  corruga¬ 
tion  of  the  waveguide  surface  can  be  introduced  lead¬ 
ing  to  scattering  of  the  propagating  fields  [21,  22].  Spa¬ 
tial  distributions  of  the  electric-field  amplitudes  of  the 
nonlinearly  interacting  fields  are  modified  this  way  and 
higher  conversion  efficiencies  of  the  nonlinear  process  can 
emerge  [23,  24].  Also  squeezed-light  generation  is  sup¬ 
ported  in  this  geometry,  as  discussed  in  [25]  where  scat¬ 
tering  of  the  second-harmonic  field  on  the  corrugation  has 
been  neglected.  A  waveguiding  geometry  with  a  periodic 
corrugation  can  also  be  conveniently  used  for  second- 
harmonic  generation  in  Cerenkov  configuration  [26].  A 
general  model  of  squeezed-light  generation  in  nonlinear 
photonic  structures  has  been  developed  in  [27]. 

In  the  pulsed  regime  and  travelling-wave  configura¬ 
tion,  squeezed-light  generation  in  the  considered  nonlin¬ 
ear  interaction  has  been  studied  with  the  help  of  phase- 
space  quasi-distributions  or  the  corresponding  Langevin 
stochastic  equations  [28].  It  has  also  been  shown  that 
there  exist  correlations  between  intensities  of  the  inter¬ 
acting  fields  in  their  transverse  planes  [29,  30].  A  local 
oscillator  in  the  form  of  an  ultrashort  optical  pulse  is 
needed  to  observe  pulsed  squeezing  experimentally.  The 
effort  to  observe  the  largest  possible  values  of  squeez¬ 
ing  has  raised  the  question  about  an  optimum  shape 
of  the  local-oscillator  field  [31,  32].  A  general  method 
how  to  solve  this  problem  is  to  use  the  Bloch-Messiah 
reduction  of  the  Green  functions  characterizing  the  lin¬ 
earized  operator  amplitude  evolution  [33].  This  method 
has  been  elaborated  in  detail  in  [34,  35]  for  degener¬ 
ate  parametric  down-conversion  in  a  BBO  crystal  and 
also  relation  between  the  Bloch-Messiah  reduction  and 
Schmidt  decomposition  of  a  two-photon  amplitude  char¬ 
acterizing  spontaneous  parametric  down-conversion  has 
been  found.  Mode  structure  of  nonclassical  states  aris¬ 
ing  from  squeezed  states  after  post-selection  using  on/off 
detectors  has  been  analyzed  in  [36]. 

In  this  report,  we  show  that  scattering  of  the  inter¬ 
acting  fields  caused  by  a  linear  periodic  corrugation  sup¬ 
ports  the  generation  of  squeezed  light  in  both  the  pump 
and  second-subharmonic  fields  provided  that  at  least  one 
of  the  interacting  fields  is  spatially  localized  as  a  conse- 
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quence  of  scattering.  We  note  that  spatial  localization 
of  optical  fields  in  this  case  is  weak  in  comparison  with 
that  occurring  in  layered  photonic  band-gap  structures 
[37-39], 

The  report  is  organized  as  follows.  In  Sec.  II,  a  quan¬ 
tum  model  of  the  nonlinear  interaction  including  both 
Heisenberg  equations  for  operator  electric-field  ampli¬ 
tudes  and  model  of  a  generalized  superposition  of  signal 
and  noise  are  presented  for  cw  pumping.  Conditions  for 
an  efficient  squeezed- light  generation  are  derived  in  Sec. 
III.  A  detailed  analysis  of  a  waveguide  made  of  LiNbC>3 
is  contained  in  Sec.  IV  considering  cw  pumping.  Gener¬ 
alization  of  the  model  to  the  case  of  multi-mode  fields  is 
given  in  Sec.  V.  Assuming  a  linear  corrugation  present 
in  the  strong  pump  field  and  weak  nonlinear  interaction 
and  relying  on  perturbation  approach,  eigenmodes  of  the 
generated  pulsed  second-subharmonic  field  are  found  in 
Sec.  VI.  Neglecting  frequency  dependence  of  the  coupling 
constants  and  inter-mode  dispersion,  the  corresponding 
equations  for  operator  amplitudes  are  solved  in  Sec.  VII. 
Numerical  analysis  of  the  principal  squeeze  variances  of 
characteristic  eigenmodes  is  contained  in  Sec.  VIII.  Con¬ 
clusions  are  drawn  in  Sec.  IX.  Appendix  A  is  devoted 
to  mode  analysis  of  an  anisotropic  planar  waveguide.  An 
optimum  mode  profile  with  respect  to  pulsed  squeezing 
is  found  in  Appendix  B. 


II.  CW  MODEL  OF  THE  NONLINEAR 
INTERACTION 

An  overall  electric-field  amplitude  E(r,f)  describing 
an  optical  field  in  the  considered  anisotropic  nonlinear 
waveguide  (shown  in  Fig.  1)  [40]  is  composed  of  two 
contributions;  pump  (or  fundamental)  electric- field  am¬ 
plitude  Ep(r,t)  at  frequency  u>  and  second-subharmonic 
electric-field  amplitude  Es(r,t)  at  frequency  ui/ 2;  i.e. 
E  =  Ep  +  Es .  We  note  that  in  case  of  second-harmonic 
generation,  the  field  with  the  amplitude  Ep(r,t)  is  called 
second-harmonic  and  that  with  the  amplitude  Es(r,f) 
is  known  as  pump.  We  keep  the  terminology  used  for 
second-subharmonic  generation  throughout  the  report. 
The  electric-field  amplitude  E  obeys  the  wave  equation 
inside  the  waveguide  with  a  nonlinear  source  term  [41] : 

o2-p  n2  p 

AE-VtV-E)-^—  =  (1) 

In  Eq.  (1)  p  denotes  vacuum  permeability,  eo  vacuum 
permittivity,  and  P„i  describes  nonlinear  polarization  of 
the  medium.  The  symbol  A  stands  for  Laplace  opera¬ 
tor,  (•)  means  a  scalar  product,  and  •  denotes  tensorial 
multiplication.  Every  spectral  component  of  the  element 
6y(r,  w)  of  the  relative  permittivity  tensor  in  the  consid¬ 
ered  waveguide  can  be  expressed  as  follows: 

dj  (x,  y,  z,  u)  =  €ij  (x,  y,  w)  [1  +  Aeij  (x,  y ,  z,  w)] .  (2) 

Small  variations  of  permittivity  e  described  by 
A eij(x,y,  z,u>)  are  induced  by  a  periodic  corrugation  of 


FIG.  1:  Four  optical  fields  interact  in  a  nonlinear  waveguide 
of  thickness  t  and  length  L  with  a  periodically-poled  (pe¬ 
riod  A„i)  x ^  susceptibility;  APf,  APb,  A„f,  and  ASb  mean 
forward-propagating  pump,  backward-propagating  pump, 
forward-propagating  second-subharmonic,  and  backward- 
propagating  second-subharmonic  electric-field  amplitudes,  re¬ 
spectively.  A  linear  corrugation  with  depth  ti  and  period  A; 
is  fabricated  on  the  waveguide  upper  surface.  Profile  of  the 
waveguide  in  the  x  —  y  plane  is  rectangular  with  width  A y 
and  depth  £;  Ay  S>  t  is  assumed.  The  waveguide  is  made  of 
LiNbC>3  with  the  optical  axis  oriented  along  the  x  axis. 


the  waveguide  surface.  These  variations  of  the  elements 
A €ij(x,y,  z,u>)  along  the  z  axis  can  be  conveniently  de¬ 
composed  into  harmonic  waves: 

OO 

A  eij(x,y,z,uj)  =  Y  eijtq(x,  y,  u)  exp 

q=—oo 

where  eq^j  are  coefficients  of  the  decomposition  and  A ; 
is  a  period  of  the  linear  corrugation.  Amplitude  of  the 
nonlinear  polarization  Pni  of  the  medium  is  determined 
using  tensor  d  of  the  second-order  nonlinear  coefficient: 

P ni(r,  t)  =  2e0d(r)  •  E(r,  t)E(r,  t).  (4) 

Taking  into  account  geometry  of  our  waveguide  elements 
djjfc(r)  of  the  nonlinear  coefficient  can  be  expressed  as 
follows: 


2tt 

iqA:Z 


,  (3) 


x,y,z)=  Y  dijk,q(x,  y)  exp 

q= — oo 

where  A n;  describes  the  period  of  a  possible  periodical 
poling  of  the  nonlinear  material. 

The  electric-field  amplitudes  of  pump  (Ep)  and  second- 
subharmonic  (Es)  monochromatic  waves  can  be  ex¬ 
pressed  in  the  form: 

Ea(a :,y,z,t)  =  i\AaF(z)ea(x,y)  exp(i/3az  -  iu>at) 

+  AaB(z)ea(x,y)  exp(-i(3az  -  iu>at) 

— H.c.] ,  a  =  p,s ,  (6) 

where  the  symbols  ep  and  es  refer  to  mode  functions  in 
the  transverse  plane  of  the  beams.  Amplitudes  APf  and 
ASf  [APb  and  ASb\  describe  forward-  [backward-]  propa¬ 
gating  pump  and  second-subharmonic  fields  and  are  such 
that  the  quantities  |APf|2,  \APb\2,  \ASf\2,  and  |ASb|2 
give  directly  the  number  of  photons  in  these  fields.  Sym¬ 
bol  /3a  means  a  propagation  constant  along  the  z  axis  in 
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mode  a  whereas  u)a  stands  for  frequency  of  this  mode. 
Symbol  H.c.  replaces  hermitian  conjugated  terms.  We 
note  that  the  considered  anisotropic  waveguide  supports 
only  TM  guided  modes;  TE  modes  are  not  guided  and 
thus  do  not  contribute  significantly  to  nonlinear  interac¬ 
tion. 

Mode  functions  ep  and  es  describe  the  transverse  pro¬ 
files  of  pump  and  second-subharmonic  fields,  respectively, 
and  fulfill  the  following  equations: 

V  x  (V  x  [ea(x,  y)  exp(±if3az  —  iu>at)]) 

U!2 

- f- €ij(x ,  y,  LOa)  ■  [ea(cc,  y)  exp (±i/3az  -  iujat )]  =  0, 

a  =  p,s.  (7) 


[42,  43] .  A  scheme  that  allows  to  decompose  interactions 
in  this  coupler  into  a  sequence  of  fictitious  interactions 
has  been  suggested  in  [44] . 

Phase  mismatches  ,  8P,  and  8niA  occurring  in  Eq.  (8) 
are  given  as  follows: 

2tt 

8 a  =  2/30  -  — ,  a  =  p,s, 

A/ 

27 r 

8ni,q  =  /3p  —  2/3s  +  q— — .  (10) 

Coefficient  q  equals  ±1  for  a  periodically  poled  nonlinear 
material,  whereas  q  =  0  for  a  material  without  period¬ 
ical  poling.  Linear  coupling  constants  Kp  and  Ks  are 
determined  along  the  expressions: 


We  note  that  mean  values  of  permittivities  (x,  y,  cua) 
are  used  for  the  determination  of  mode  functions.  The 
mode  functions  ep  and  es  are  normalized  to  describe  one 
photon  in  a  mode  inside  the  waveguide  (see  Appendix 

A). 

Substitution  of  Eqs.  (2-6)  into  Eq.  (1)  assuming 


d2Aa 


dz2 


< 


Pa 


OAa 


dz 


for  a  =  p,s  and  b  =  F,  B  (analog 


of  the  slowly-varying  envelope  approximation  to  spatial 
evolution)  results  in  the  following  equations  for  ampli¬ 
tudes  APf  ,  APb  ,  ASf  ,  and  ASB : 


dASF 

dz 

dASB 

dz 

dAPF 

dz 

dAPB 

dz 


iKs  exp (-i8sz)ASB 
T  8:Kni,q  &xp(%8ni,qz}  ApF  ASF , 

-iK*  exp(i(5sz)ASF 

8Kni.q  exp(  idni^qz')  APb  ASB , 
iKp  exp(—i8pz)APB 
-  2Kni,qexP(-idm,qz)A2SF, 

— iK *  exp(iSpz)APF 

+  2Kni,q  exp (i8ni,qz)A2SB .  (8) 


Because  the  waveguide  is  made  of  an  anisotropic  ma¬ 
terial,  also  the  following  relation  has  to  hold  in  order  to 
derive  correctly  equations  in  Eq.  (8)  (for  details,  see  [41]): 


d^p^[ea{x,y)]z  exp(±ipaz)  z 
rIA  M 

+  — ([ea(x,  y)}z  exp [±ipaz)) 

+  dA°,b ^  (V  •  [ea(x,  y)  exp(±ipaz)])  z  «  0, 
dz 

a  =  p,s,  b  =  F,B,  (9) 


Kn  = 


wa  I  dxdye-i{x,y,LOa)  ■  e*a(x,y)ea(x,y) 


2  c2pa 
a  =  p,s. 


f  dxdy  \ea(x,y)\2 


(11) 


Similarly,  the  following  expression  can  be  found  for  non¬ 
linear  coupling  constants  Knitq  for  q  =  0,  ±1: 


2i(x2Jdxdydq(x,y)  ■  ep(x,y)e*s(x,y)e*s(x,y) 
c2Ps  J  dxdy\es(x,y)\2 

2iu2  f  dxdy  dq(x,y)  ■  ep(x,y)e*s(x,y)e*s(x,y) 
c2Pp 


f  dxdy  \ep(x,y)\2 


(12) 


The  last  approximate  equality  in  Eq.  (12)  is  valid  pro¬ 
vided  that  cos  IPs  ~  ojp/Pp  and  due  to  the  normalization 
of  the  mode  functions  es  and  ep.  This  approximation  as¬ 
sures,  that  only  one  nonlinear  coupling  constant  occurs 
in  Eqs.  (8)  which  is  important  in  quantum  description. 

Expressions  for  linear  ( Ks ,  Kp)  and  nonlinear  (Kni  0, 
Kni,i)  coupling  constants  appropriate  for  the  considered 
waveguide  and  derived  from  Eqs.  (11)  and  (12)  can  be 
found  in  Appendix  A  in  Eqs.  (A10-A12). 

Quantum  model  of  the  nonlinear  interaction  in  the  con¬ 
sidered  waveguide  can  be  formulated  changing  the  clas¬ 
sical  envelope  electric-field  amplitudes  APf,  APb,  ASf, 
and  ASb  occurring  in  Eq.  (6)  into  operators  denoted  as 
APf,  Apb,  ASf,  and  ASB,  respectively.  A  quantum  ana¬ 
log  of  the  classical  equations  written  in  Eq.  (8)  can  then 
be  derived  from  the  Heisenberg  equations  (for  details,  see 
[45,  46]; 


dX 

dz 


i 

n  l 


G,X 


(13) 


considering  the  following  momentum  operator  G: 


G  = 


hKsexp(-idsz)Al  A 


Sf  sb 


and  z  stands  for  a  unit  vector  along  the  z  axis.  We  note 
that  Eqs.  (8)  describe  also  a  nonlinear  coupler  composed 
of  two  waveguides  made  of  y*-2-*  media  whose  modes  in¬ 
teract  through  evanescent  waves.  Nonclassical  properties 
of  light  propagating  in  this  coupler  have  been  studied  in 


+hKpexp(—i8pz)A^pAPB  +H.c. 

2ihKni,q  exp(i8nitqz)A\pAPF 
+2ihK„iiq  exp(  i8ni  qz)A^SB  -Apb  H.c.j  ,  (14) 
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where  (7  =  0org  =  ±l.  In  Eq.  (13),  h  means  the  reduced 
Planck  constant,  X  stands  for  an  arbitrary  operator  and 
symbol  [,  ]  denotes  commutator. 

The  operator  quantum  equations  analogous  to  those 
written  in  Eq.  (8)  can  be  solved  using  the  method  of 
a  small  operator  correction  (denoted  as  6 A)  to  a  mean 
value  (denoted  as  A)  in  which  an  electric-field  opera¬ 
tor  amplitude  A  is  decomposed  as  A  =  A  +  6 A.  This 
method  provides  a  set  of  classical  nonlinear  equations  for 
the  mean  values  A  that  coincides  with  the  set  given  in 
Eq.  (8) .  The  operator  electric-field  amplitude  corrections 
5 A  fulfill  the  following  linear  operator  equations: 

apb6A*sb  +  a*sb6Apb]  , 

5ASf  , 

1sJASb.  (15) 

Functions  JCS,  Kpi  and  Kni  q  introduced  in  Eqs.  (15)  are 
defined  as: 

Ka  =  iKaexp(-iSaz),  a=s,p, 

ICF,q  =  4Kni,qexp(i5ni,qz), 

ICB,q  =  4:Kni,qexp(-iSnitqz).  (16) 


d6Aa 


dz 

d8ASB 

dz 

dSAPF 

dz 

dSAVo 


dz 


K-s&Asb  +  K,F,q  | 
X*s5ASf  -  KB,q 
ICpSAPB  —  K,*p qj 


K.pSApB  +  X*B  q. 


We  note  that  any  solution  of  the  nonlinear  equations  in 
Eqs.  (8)  obeys  the  following  relation  useful  in  a  numerical 
computation: 


^  {\ASF  |2  +  2| I2  -  |^SB  |2  -  2|APb  I2)  =  0.  (19) 


The  solution  of  the  system  of  linear  operator  equations 
in  Eqs.  (15)  for  the  operator  electric-field  amplitude  cor¬ 
rections  5 A  can  be  found  numerically  and  put  into  the 
following  matrix  form: 


f  Mb, out  \  _  (  UfF  Ufb  \  f  Mgin  \ 
V  Mb, in  /  \Ubf  Ubb  )  \  Mb, out/ 

where 


(20) 


Mb, in  = 


/SASF(  0)\ 

mUo) 

^f(O) 

\sAIf(o)J 


8Ab,  in 


fSA3B(L)\ 
sA\B(L) 
&Apb  (l) 
\sAIb(l)J 


6Af , 


out  — 


fSASF(L)\ 
sA\F(L) 
fiAPF  ( L ) 
\6AlF(L)J 


^-^BjOUt 


(6Asb(  0)\ 

mU°) 

fiAPB  (o) 
\sAIb(0)J 


(21) 


and  L  means  length  of  the  waveguide.  Matrices  Uff, 
Ufb,  Ubf,  and  UBb  are  determined  using  numerical  so¬ 
lution  of  Eqs.  (15). 

The  following  input-output  relations  among  the  oper¬ 
ator  amplitude  corrections  8 A, 


Solution  of  the  classical  nonlinear  equations  written  in 
Eqs.  (8)  can  only  be  reached  numerically  using,  e.g.,  a 
finite  difference  method  called  BVP  [47].  This  method 
requires  an  initial  guess  of  the  solution  that  can  be  con¬ 
veniently  obtained  when  the  nonlinear  terms  in  Eqs.  (8) 
are  omitted.  Then,  the  initial  solution  can  be  written  as 
follows: 


Aft 


A<0) 

aB 


=  exp  I  —i 


Ba  cos(Aaz)  +  Ba  sin(Aa,z) 


=  exp  l  i 


Sn.Z 


Ba  \  -7F7F-  cos(Aa2)  +  sin(A0z) 
\  ^-K-a  -K-a 

+Ba  \-jA^sm(Aaz)  -i^cos(Aaz) 

a  =  s,p, 


(17) 


and 


2 

a  |  > 


a  =  5,p. 


(18) 


In  Eqs.  (17),  constants  Bp,  Bp ,  Bs,  and  Bs  are  set  ac¬ 
cording  to  the  boundary  conditions  at  both  sides  of  the 
waveguide. 


f  ^^BjOut  \ 

V  ^^4s,out  / 


(  Uff  —  UfbUbbUBf 
\  ~UbbUbf 

x  (  SAf^  ) 

V  Ms, in  / 


UfbUbb 

Ubb 


(22) 

(23) 


are  found  solving  Eqs.  (21)  with  respect  to  vectors 
Mb, out  and  Mb, out-  The  output  operator  electric-field 
amplitude  corrections  contained  in  vectors  Mb, out  and 
Mb, out  obey  boson  commutation  relations  provided  that 
the  input  operator  electric-field  amplitude  corrections 
given  in  vectors  Mb, in  and  Mb, in  are  ruled  by  boson 
commutation  relations.  We  note  that  also  certain  com¬ 
mutation  relations  among  the  operator  electric-field  am¬ 
plitude  corrections  in  vectors  Mb, out  and  Mb, in  can  be 
derived  (for  details,  see  [46]). 

We  restrict  our  considerations  to  states  of  optical  fields 
that  can  be  described  using  the  generalized  superposition 
of  signal  and  noise  [2].  Thus  coherent  states,  squeezed 
states  as  well  as  noise  can  be  considered.  Parameters 
Bj ,  Cj,  Djk,  and  Dj k  defined  below  are  sufficient  for  the 
description  of  any  state  of  a  two-mode  optical  field  in  this 
approximation  [20]: 

Bj  =  (AijAi,-), 
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Cj  =  ((A  A3f), 

Djk  =  (AAjAAk),  j  ^  k, 

Djk  =  -(AitAifc),  jVfc;  (24) 

AAj  =  Aj  —  ( Aj );  symbol  (  )  denotes  the  quantum 

statistical  mean  value.  Expressions  for  the  coefficients 
Bj ,  Cj,  Djk,  and  Djk  appropriate  for  outgoing  fields  can 
be  derived  [46]  using  elements  of  matrix  U  given  in  Eq. 
(23)  and  incident  values  of  coefficients  BjjTIA  and  C:h\n  ^ 
related  to  anti-normal  ordering  of  field  operators  (for  de¬ 
tails,  see  [20]): 

Djjn.A  =  cosh  ( T j )  “h  IT'chJ, 

Cj,. in, a  =  ^  exp(ii?j)  sinh(2rj).  (25) 

In  Eq.  (25),  Tj  stands  for  a  squeeze  parameter  of  the 
incident  j-th  field,  dj  means  a  squeeze  phase,  and  nch,j 
stands  for  a  mean  number  of  incident  chaotic  photons. 
Coefficients  Djkjn,A  and  Djkjn,A  for  an  incident  field  are 
considered  to  be  zero,  i.e.  the  incident  fields  are  assumed 
to  be  statistically  independent. 

The  maximum  attainable  value  of  squeezing  of  electric- 
field  amplitude  fluctuations  is  given  by  the  value  of  a 
principal  squeeze  variance  A  [2].  Both  single-mode  prin¬ 
cipal  squeeze  variances  A  j  and  compound-mode  principal 
squeeze  variances  A ij  (characterizing  an  overall  field  com¬ 
posed  of  two  other  fields)  can  be  determined  in  terms  of 
coefficients  Bj,  Cj,  Djk,  and  Djk  given  in  Eq.  (24)  (for 
details,  see,  [2,  20]): 

Xj  =  l  +  2[Bj-\Cj\\,  _  (26) 

A jk  =  2  [l  +  Bj  +  Bk  —  2 ^R(Djk) 

-\Cj  +  Ck  +  2Djk\}-  (27) 

symbol  3?  means  the  real  part  of  an  expression.  Values 
of  the  principal  squeeze  variance  A  j  ( \jk )  less  than  one 
(two)  indicate  squeezing  in  a  single-mode  (compound¬ 
mode)  case. 

III.  SUITABLE  CONDITIONS  FOR 
SQUEEZED-LIGHT  GENERATION 

It  occurs  that  there  exist  two  conditions  for  an  ef¬ 
ficient  squeezed-light  generation.  The  first  condition 
comes  from  the  requirement  that  the  nonlinear  interac¬ 
tion  should  be  phase-matched  along  the  whole  waveguide, 
whereas  the  second  one  gives  optimum  conditions  for  the 
localization  of  interacting  optical  fields. 

A.  Overall  phase-matching  of  the  nonlinearly 
interacting  fields 

Conditions  for  an  optimum  phase-matching  of  the  in¬ 
teracting  fields  can  be  revealed,  when  we  write  the  dif¬ 
ferential  equation  for  the  number  Na  of  photons  in  field 


a;  Na  =  A*Aa.  Using  Eqs.  (8)  we  arrive  at  the  following 
differential  equations: 

-29  {ICS  exp(-iSsz)A*SF  ASb  } 

+  8$  {Knkq  exp(i6ni,qz)A*s2FAPF}  , 

-29  { Ks  exp(-iSsz)A*SF  ASB  } 

-  83?  {KnUq  exp  (~i6nhqz)A*s2BAPB  }  , 

-29  {Kp  exp (-iSpz)A*pFAPB  } 

-  43?  {KnUq  exp(iSnlqz)A*s2FAPF  }  , 

-29  {Kp  exp (-iSpz)ApFAPB  } 

+  m{Knl:qexV(-iSnl:qz)AfBAPB}-,  (28) 

symbol  9  denotes  the  imaginary  part  of  an  expression. 
We  note  that  the  first  terms  on  the  right-hand  side  of 
the  first  and  the  second  (as  well  as  the  third  and  the 
fourth)  equations  in  Eqs.  (28)  have  the  same  sign  because 
of  counter-propagation  of  the  fields  [see  also  Eq.  (19)]. 
The  nonlinear  interaction  described  by  the  second  terms 
on  the  right-hand  sides  of  Eqs.  (28)  is  weak  and  so  we  can 
judge  the  contribution  of  these  terms  using  a  perturba¬ 
tion  approach.  In  the  first  step  we  neglect  the  nonlinear 
terms  in  Eqs.  (8)  and  solve  Eqs.  (8)  for  field  amplitudes 
Agp,  A^,  A^Pp,  and  A P°J.  Then  we  insert  this  solution 
into  the  nonlinear  terms  in  Eqs.  (28)  and  find  this  way 
optimum  conditions  that  maximize  contributions  of  these 
terms.  The  solution  for  amplitudes  AiB ,  A^B,  A{0J,  and 
Ap°J  coincides  with  that  written  in  Eqs.  (17)  as  an  initial 
guess  for  numerical  solution  and  we  rewrite  it  into  the 
following  suitable  form: 

•W)  =  exp 

X  [Kf  exp(?'Aaz)  +  B~f  exp(— iA0z)]  , 

y,a  B a  iDa 

VaF  -  ^  ’ 

_  Ba  +  iBa 
VaF  ~  ^  ’ 

Aa “2  (*)  =  eXp(*^) 

x  [Kb  exp(?'Aaz)  +  B~b  exp(— zAaz)]  , 
B:B  =  ~\+K2aAa(Ba-iBa), 

BaB  =  +  iBa),a  =  p,  s. 

(29) 

Provided  that  a  linear  corrugation  is  missing  in  field  a 
the  solution  for  amplitudes  AaF  and  A^B  can  be  obtained 
from  the  expressions  in  Eqs.  (29)  using  a  sequence  of  two 
limits;  Sa  — >  0,  Ka  — >  0: 

=  Ba, 


dNSF 

dz 

dNSB 

dz 

dNPF 

dz 

dNPB 

dz 
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42(*)  =  Ba.  (30) 

Interaction  between  the  forward-propagating  fields  is 
described  in  Eqs.  (28)  in  our  perturbation  approach  by 
the  term 

3?  exp  (i6ni,qz)A(°J*2A(°J  J  (31) 

that,  after  substituting  the  expressions  for  amplitudes 
and  Ap^  from  Eqs.  (29),  splits  into  the  following 
eight  terms: 

»{  KnliVB?/B?/BfF 

x  exp (i[Sni,q  +  Ss  -  Sp/2  =f  As  =f  As  ±  Ap]z)  j.(32) 


The  nonlinear  interaction  is  efficient  under  the  condition 
that  one  of  these  terms  does  not  oscillate  along  the  z 
axis.  This  gives  us  eight  possible  conditions  that  com¬ 
bine  nonlinear  phase-mismatch  Sni^q  and  parameters  of 
the  corrugation  S3,  6p,  Ks ,  and  Kp: 

Sni,q  +  Ss-  Sp/2  =f  As  =f  As  ±  Ap  =  0.  (33) 

It  depends  on  a  given  waveguide  and  initial  conditions 
which  out  of  these  eight  conditions  leads  to  an  efficient 
nonlinear  interaction.  We  note  that  the  conditions  in 
Eq.  (33)  are  valid  also  for  the  nonlinear  interaction  be¬ 
tween  the  backward-propagating  fields  characterized  by 

the  term  3?  ^Kni,q  exp (-iSni,qz)Ai°J  AP°J  j  in  Eqs.  (28). 

We  consider  two  special  cases  in  which  a  linear  cor¬ 
rugation  is  present  either  in  the  pump  or  the  second- 
subharmonic  field.  Assuming  the  linear  corrugation  in 
the  pump  field  the  conditions  in  Eq.  (33)  get  the  form: 


6ni,g  -  Sp/2  ±  Ap  =  0.  (34) 

Sign  -  (+)  is  suitable  for  SniA  >  0  (■ 8nitq  <  0)  when  we 
solve  Eq.  (34)  for  Sp: 


Sp  Snl  q  T 


\KP\\ 

Snl,q 


(35) 


i.e.  Snit q  and  Sp  have  the  same  sign.  The  expression  in 
Eq.  (35)  then  determines  the  period  A;  of  linear  corruga¬ 
tion: 

7T 

A/  =  PP-Snl,q/2~\Kp\y(25nl,q)-  (36) 

On  the  other  hand,  the  conditions 

Sni,q  +  Sa/2  =f  As  =F  As  =  0  (37) 

are  suitable  for  the  linear  corrugation  in  the  second- 
subharmonic  field.  When  Sniiq  >  0  (dn;,9  <  0)  signs  - 
(+)  in  Eq.  (37)  are  appropriate  and  we  have: 

Sni,Q  2|  Ks\2 


6„  =  — 


(38) 


2  5nhq 

i.e.  8ni:q  and  Ss  have  opposite  signs.  The  period  A /  of 
linear  corrugation  is  then  given  as 

tt 


A  ,= 


Ps  +  Snl,q/ 4  + 


\KP\2/Snltq 


(39) 


B.  Localization  of  the  interacting  fields 

If  we  have  considered  a  waveguide  with  the  linear  cor¬ 
rugation  so  deep  that  it  touches  the  bottom  of  the  waveg¬ 
uide,  we  would  have  a  layered  structure  with  band-gaps 
and  transmission  peaks.  If  a  field  frequency  lies  in  a 
transmission  peak,  the  field  is  also  well  localized  inside 
the  structure  [37].  The  considered  waveguide  with  a  weak 
linear  corrugation  behaves  qualitatively  in  the  same  way 
[23] .  The  areas  of  transparency  of  the  waveguide  are  con¬ 
venient  for  the  nonlinear  interaction  because  of  the  field 
localization  (though  very  weak  for  a  weak  corrugation). 
A  transmission  peak  of  the  waveguide  can  be  found  from 
the  condition  that  a  backward-propagating  field  is  zero 
both  at  the  beginning  and  at  the  end  of  the  waveguide, 
because  it  is  not  initially  seeded.  These  requirements  are 
fulfilled  by  the  solution  in  Eqs.  (17)  under  the  following 
conditions: 

mn 

A  a  =  ~j~,  m=  1,2,...; 

AX 

Ba  =  45(0);  Ba=^-Ba,  a  =  p,s.  (40) 

Natural  number  m  counts  areas  of  transmission. 

The  linear  phase  mismatch  Sa  is  determined  from  the 
first  equation  in  Eqs.  (40)  in  the  following  form 


Sa  =  ±2 \J {jY~y  +  \Ka\2,  a=p,s,  (41) 

and  the  corresponding  period  A i  of  linear  corrugation  is 
given  as: 


Pa  ±  \J{ nrni/L )2  +  | Ka \ 2 


(42) 


The  conditions  in  Eqs.  (35)  and  (41)  for  a  linear  cor¬ 
rugation  in  the  pump  field  can  be  combined  together  to 
provide  the  following  formula  for  the  coupling  constant 

Ap: 


(43) 

The  linear  phase  mismatch  Sp  is  then  determined  along 
the  formula  in  Eq.  (41)  such  that  the  signs  of  5ni,q  and 
Sp  are  the  same. 

Similarly  the  conditions  for  a  linear  corrugation  in  the 
second-subharmonic  field  written  in  Eqs.  (38)  and  (41) 
and  considered  together  lead  to  an  expression  for  the  cou¬ 
pling  constant  Ks : 


\Ks\ 2 


(44) 


The  sign  of  linear  phase  mismatch  Ss  determined  from 
Eq.  (41)  is  opposite  to  that  of  Sn/,g- 
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IV.  SQUEEZED-LIGHT  GENERATION  - 
NUMERICAL  ANALYSIS  FOR  CW  PUMPING 


Discussion  of  squeezed-light  generation  is  decomposed 
into  three  parts.  In  the  first  one,  we  pay  a  detailed  atten¬ 
tion  to  attainable  characteristics  of  the  considered  waveg¬ 
uide.  The  second  part  is  devoted  to  second-subharmonic 
generation,  i.e.  an  incident  strong  pump  field  is  assumed. 
In  the  third  part,  a  strong  incident  second-subharmonic 
field  is  assumed,  i.e.  second-harmonic  generation  is  stud¬ 
ied. 


A.  Characteristic  parameters  of  the  waveguide 

We  consider  a  waveguide  made  of  LiNbC>3  with  length 
L  =  lx  1CT3  m  and  width  Ay  =  1  x  10-5  m  pumped 
at  the  wavelength  Xp  =  0.534  x  10-6  m  (As  =  1.064  x 
10-6  m).  More  details  are  contained  in  Appendix  A. 
We  require  a  single-mode  operation  at  both  the  pump- 
and  second-subharmonic-field  frequencies  that  can  be 
achieved  for  the  thickness  t  of  the  waveguide  in  the  range 
t  e  (0.403,0.544)  x  10-6  m  [24].  We  also  neglect  losses 
in  the  waveguide  that  may  be  caused  both  by  absorption 
inside  the  waveguide  and  scattering  of  light  that  does  not 
propagate  into  guided  modes. 

Values  of  the  natural  nonlinear  phase  mismatch  8nit o 
are  high  (around  1.7  x  106  nU1)  in  the  region  of  single¬ 
mode  operation  and  depend  on  the  thickness  t  of  the 
waveguide  (see  Fig.  2a).  Values  of  the  nonlinear  cou¬ 
pling  constant  A'„/,o  depicted  in  Fig.  2b  increase  with 
the  increasing  values  of  the  thickness  t  because  of  the 
thicker  the  waveguide  the  more  the  modes  are  localized 
inside  the  waveguide  and  so  the  greater  the  overlap  of  the 
nonlinearly  interacting  field  amplitudes.  Attainable  val¬ 
ues  of  linear  coupling  constants  Kp  and  Ks  for  different 
values  of  thickness  t  and  depth  of  corrugation  ti  can  be 
obtained  from  Fig.  3.  Values  of  the  coupling  constants 
Kp  and  Ks  for  a  given  thickness  t  increase  monotonously 
with  the  increasing  values  of  the  depth  of  corrugation  ti. 
For  small  values  of  the  thickness  t,  values  of  the  coupling 
constant  Ks  are  small  because  the  waveguide  is  thin  for 
the  second-subharmonic  field  and  so  a  considerable  part 
of  the  field  is  outside  the  waveguide  and  cannot  be  scat¬ 
tered  by  the  corrugation. 

Amplitudes  of  the  dimensionless  incident  strong 
electric-field  amplitudes  Aah  are  determined  from  the  in¬ 
cident  power  Pah  along  the  relation: 


\Aab\  =  a  =  P,s,  b=F,B.  (45) 

On  the  other  hand,  power  P° kut  of  an  outgoing  field  is 
given  as  follows: 
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FIG.  2:  Natural  nonlinear  phase  mismatch  5ni, o  (a)  and  ab¬ 
solute  value  of  the  nonlinear  coupling  constant  Kni, o  (b)  as 
they  depend  on  the  thickness  t  of  the  waveguide  are  shown  in 
the  region  of  single-mode  operation;  constant  Kni,o  is  deter¬ 
mined  for  amplitudes  that  correspond  to  one  pump  and  one 
second-subharmonic  photon  inside  the  waveguide. 

=  -KfNab,  a  =  p,s,  b  =  F,B ;  (46) 

Pa-LJ 

Nab  denotes  the  number  of  photons  leaving  the  waveg¬ 
uide. 

New  dimensionless  parameters  are  convenient  for  the 
discussion  of  behavior  of  the  waveguide; 

-  z/L, 

A[  =  AJ/L,A^  =  A  niJL, 

Pra  =  Lpa,  Sra  =  LSa,  K'a  =  LKa,  a=p,s, 

5rnhq  =  L5nl,q,Krnl:q  =  LKnltq.  (47) 

Applying  these  parameters  the  waveguide  extends  from 
zr  =  0  to  zr  =  1.  The  dimensionless  parameters  en¬ 
able  to  understand  the  behavior  of  the  waveguide  as  it 
depends  on  the  length  L  using  graphs  and  discussion  bel¬ 
low. 


B.  Second-subharmonic  generation 

As  a  reference  for  the  efficiency  of  squeezed-light  gen¬ 
eration  we  consider  the  waveguide  with  periodical  poling 
and  assume  that  it  is  pumped  by  the  power  of  2  W.  The 
nonlinear  interaction  is  perfectly  phase  matched  for  the 
period  of  poling  h.rnl  ss  3.552  x  1CU3  where  we  have  for 
the  principal  squeeze  variances  XSp  s=s  0.45  and  XPF  ~  1 
(see  Fig.  4a).  The  more  distant  the  value  of  Krnl  from  the 
above-mentioned  optimum  value  is,  the  larger  the  non¬ 
linear  phase  mismatch  and  the  larger  the  value  of  the 
principal  squeeze  variance  XSF. 


a) 


0.0380 

0.0355 

0.0329 

0.0304 

0.0279 

0.0253 

0.0228 

0.0203 

0.0177 

0.0152 

0.0127 

0.0101 

0.0076 

0.0051 

0.0025 

0 


0.0070 
0.0065 
0.0060 
0.0055 
0.0050 
0.0045 
0.0040 
0.0035 
0.0030 
0.0025 
0.0020 
0.0015 
0.0010 
0.0005  l 
0 

t  (x10'6m) 


FIG.  3:  Contour  plot  of  the  absolute  value  of  pump-field 
[second-subharmonic-field]  linear  coupling  constant  Kp  [Ks\ 
as  it  depends  on  thickness  t  of  the  waveguide  and  depth  ti 
of  linear  corrugation  is  shown  in  (a)  [(b)]  in  the  region  of 
single-mode  operation. 


Considering  a  periodic  corrugation  in  the  pump  field, 
better  values  of  squeezing  in  the  second-subharmonic 
field  can  be  reached.  The  reason  is  that  the  corruga¬ 
tion  causes  scattering  of  the  pump  field  that  leads  to 
inhomogeneity  of  the  field  along  the  2  axis.  The  nonlin¬ 
ear  process  then  profits  from  field  localization.  However, 
nonzero  values  of  the  nonlinear  phase  mismatch  Snit i  are 
important  to  exploit  this  localization  because  they  have 
to  compensate  periodic  spatial  oscillations  caused  by  the 
corrugation  (see  [48]).  A  perfect  phase  matching  of  all 
the  processes  occurring  in  the  waveguide  can  be  reached 
this  way. 

A  typical  dependence  of  the  principal  squeeze  vari¬ 
ance  XSF  as  well  as  the  number  NSf  of  photons  leav¬ 
ing  the  waveguide  for  the  forward-propagating  second- 
subharmonic  field  for  attainable  values  of  parameters  of 
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FIG.  4:  Principal  squeeze  variances  \SF  (solid  curve)  and 
A PF  (solid  curve  with  *)  as  functions  of  the  period  A ^  of  pol¬ 
ing  for  second-subharmonic  generation  ( PVF  =  2  W,  Psp  = 
1  x  10-10  W  -  a  negligible  seeding  to  substitute  spontaneous 
process  in  classical  equations)  (a)  and  second-harmonic  gener¬ 
ation  ( PPF  =  1  x  1CP10  W,  PSp  =  2  W)  (b);  to  understand  the 
behavior  with  respect  to  the  length  L,  the  period  A should 
be  decomposed  as  A°ni  +AA„i,  then  5^  q  «  — 2-7rqAA[u/A°i  , 
A  =  2nq/(2(3a  -  &),  AArnl  =  LAAni ;  (  =  5x  1(T7  m, 
I<p  =  0,  Krs  =  0,  arg(App)  =  0,  arg(Asp)  =  0,  APB  =  0, 
Asb  =  0;  A°nl  «  3.552  x  1CT6  nr. 


the  corrugation  is  shown  in  Fig.  5  assuming  a  fixed  value 
of  the  nonlinear  phase  mismatch  S'nl  1  equal  to  -10.82 
(Arnl  =  3.53  x  10-3),  i.e.  it  corresponds  roughly  to  the 
first  local  minimum  of  XSF  in  the  curve  in  Fig.  4a.  We 
can  clearly  see  that  an  efficient  nonlinear  interaction  oc¬ 
curs  in  strips  that  correspond  to  transmission  peaks;  the 
larger  the  number  m  of  a  transmission  peak  [see  Eq.  (40)] 
the  weaker  the  effective  nonlinear  interaction.  The  neces¬ 
sity  to  fulfill  also  the  condition  for  perfect  overall  phase 
matching  given  in  Eq.  (35)  is  evident.  The  principal 
squeeze  variance  XSF  reaches  values  around  0.2  inside 
the  strips  around  the  first  several  transmission  peaks. 
The  larger  the  number  to  of  a  transmission  peak,  the 
greater  values  of  linear  coupling  constant  K'j  and  linear 
phase-mismatch  Srp  have  to  be  used  to  reach  high  levels 
of  squeezing.  Up  to  several  forward-propagating  pho¬ 
tons  can  be  present  inside  the  waveguide  (see  Fig.  5b) 
at  a  given  time  instant.  This  means  that  the  power  of 
the  outgoing  field  is  of  the  order  of  10-8  W  (energy  of 
one  second-subharmonic  photon  inside  the  waveguide  of 
thickness  t  =  5  x  10“ 7  m  and  length  L  =  lx  1CU3  m  cor¬ 
responds  to  the  output  power  of  2.58  x  10-8  W).  Only 
the  first  and  the  second  transmission  peaks  can  give  rea¬ 
sonable  values  of  the  power  of  the  outgoing  field. 

Role  of  the  value  of  nonlinear  phase  mismatch  5 ^  1  to 
squeezed-light  generation  is  revealed  in  Fig.  6,  where  the 
principal  squeeze  variance  XSF  as  a  function  of  the  non- 
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FIG.  5:  Contour  plots  of  principal  squeeze  variance  ASF  (a) 
and  number  Nsp  of  photons  leaving  the  waveguide  (b)  for 
the  forward-propagating  second-subharnronic  field  as  they  de¬ 
pend  on  parameters  \Kp\  and  5p  of  the  corrugation  in  the 
pump  field;  more-less  equidistant  curves  in  the  graphs  indi¬ 
cate  positions  of  transmission  peaks  [Eq.  (40)] ,  the  last  curve 
going  up  gives  the  condition  in  Eq.  (35);  Arnl  =  3.53  x  10-3 
(%,i  =  -10.82),  PPF  =  2  W,  PSF  =  1  x  to"10  W,  arg(it'p)  = 
7t/2;  values  of  the  other  parameters  are  the  same  as  in  Fig.  4. 
Triangles  in  lower  right  corners  with  formally  zero  values  in 
both  graphs  lie  in  the  region  with  an  exponential  behavior  of 
classical  amplitudes  that  is  not  suitable  for  nonclassical-light 
generation. 


linear  phase  mismatch  1  under  the  conditions  given  in 
Eqs.  (41)  and  (43)  is  drawn  for  the  first  two  transmission 
peaks  (to  =  1,2).  The  larger  the  value  of  5^  1  the  better 
the  values  of  the  principal  squeeze  variance  ASf  ,  but  also 
the  greater  the  values  of  linear  coupling  constant  Kp  \ 
and  linear  phase  mismatch  Sp  (as  shown  in  Figs  6b  and 
6c).  The  smallest  values  of  Kp  |  and  5p  are  obtained  when 
the  generation  occurs  around  the  first  transmission  peak. 
From  practical  point  of  view,  the  depth  of  linear  corru¬ 
gation  is  limited  by  technological  reasons  (also  validity  of 


the  model  has  to  be  judged  for  greater  values  of  the  linear 
coupling  constant  Kp )  and  this  means  that  only  smaller 
values  of  the  linear  coupling  constant  I\v  are  available 
(see  Fig.  3).  In  order  to  obtain  the  best  possible  value 
of  the  principal  squeeze  variance  ASf  the  curves  in  Fig.  6 
suggest  to  work  in  the  first  transmission  peak  and  use 
the  greatest  available  value  of  Kp  (the  deepest  possible 
corrugation).  From  the  curve  in  Fig.  6b  the  appropriate 
value  of  the  nonlinear  phase  mismatch  1  is  determined 
and  then  the  curve  in  Fig.  6c  gives  a  suitable  value  of  the 
linear  phase  mismatch  Sp  (and  consequently  the  period 
A i  of  linear  corrugation).  Keeping  the  incident  pump- 
field  power  fixed,  the  depth  f;  of  the  corrugation  limits 
the  achievable  values  of  the  principal  squeeze  variance 
ASf;  the  deeper  the  corrugation,  the  better  the  squeez¬ 
ing.  We  can  roughly  say  that  the  appropriate  values  of 
linear  coupling  constant  Kp  and  linear  phase  mismatch 
Sp  are  close  (or  higher)  to  the  value  of  nonlinear  phase 
mismatch  1  to  have  an  optimum  phase  matching  (see 
Figs.  6b, c). 

The  curves  in  Fig.  6a  also  reveal  the  role  of  length  L  of 
the  waveguide  in  squeezed-light  generation.  Because  the 
nonlinear  phase  mismatch  <5ra/ i  does  not  depend  on  the 
length  L ,  the  dimensionless  nonlinear  phase  mismatch 
<5 ^  1  depicted  in  Fig.  6a  is  linearly  proportional  to  the 
length  L  and  so  the  larger  the  value  of  length  L  the  lower 
the  value  of  the  principal  squeeze  variance  XSF .  The  anal¬ 
ysis  of  expressions  in  Eqs.  (41)  and  (43)  shows  that  in  the 
limit  of  large  length  L ,  \Kp\  — >  Sn^±  and  Sp  — >  2Snit\.  We 
note  that  the  distance  of  the  adjacent  transmission  peaks 
in  the  plane  spanned  by  variables  Kp  and  Sp  behaves  as 
1/L  (see,  e.g.,  graphs  in  Fig.  5). 

In  case  of  the  considered  LiNbC>3  waveguide,  values  of 
the  natural  nonlinear  phase  mismatch  5nito  for  the  con¬ 
sidered  values  of  thickness  t  as  given  in  Fig.  2a  are  so 
high  that  in  order  to  take  advantage  of  scattering  on  the 
corrugation,  an  additional  periodical  poling  has  to  be  in¬ 
troduced  to  partially  compensate  for  this  mismatch  [25] . 

Usefulness  of  the  corrugation  can  be  demonstrated 
when  we  compare  the  following  three  configura¬ 
tions:  perfectly  quasi-phase  matched  waveguide,  non- 
perfectly  quasi-phase  matched  waveguide,  and  finally 
non-perfectly  quasi-phase  matched  waveguide  with  a 
suitable  corrugation  that  compensates  for  phase  mis¬ 
match.  The  waveguide  with  a  corrugation  gives  better 
values  of  the  principal  squeeze  variance  XSF  and  also 
considerably  greater  values  of  the  number  NSp  of  pho¬ 
tons  leaving  the  waveguide  in  comparison  with  the  per¬ 
fectly  quasi-phase-matched  waveguide,  as  documented  in 
Fig.  7.  Also  an  improvement  caused  by  an  introduction  of 
the  corrugation  into  a  non-perfectly  quasi-phase-matched 
waveguide  is  worth  mentioning. 

A  linear  corrugation  can  be  alternatively  introduced 
into  the  second-subharmonic  field.  Despite  the  fact  that 
the  second-subharmonic  field  remains  at  single  photon 
level,  results  obtained  for  principal  squeeze  variances 
A  and  numbers  N  of  photons  are  comparable  to  those 
achieved  with  a  corrugation  in  the  pump  field.  Suitable 
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FIG.  6:  Principal  squeeze  variance  XSF  (a),  absolute  value  of 
the  linear  coupling  constant  Kp  given  by  Eq.  (43)  (b),  and 
linear  phase  mismatch  5p  determined  along  Eq.  (41)  (c)  as 
functions  of  the  nonlinear  phase  mismatch  p  first  (m=l) 
and  second  (m=2)  transmission  peaks  as  well  as  signs  +  and 
-  in  Eq.  (43)  are  considered  (m  =  1+:  solid  curve  with  *, 
m  =  1  —  :  solid  curve  with  o,  m  =  2+:  solid  curve  with  o, 
m  =  2—:  solid  curve  with  A);  PPF  =  2  W,  Pap  =  lx  10_1°  W, 
arg (Kp)  =  7r/2;  values  of  the  other  parameters  are  the  same 
as  in  Fig.  4. 


conditions  are  given  in  Eqs.  (38)  and  (40)  in  this  case. 


C.  Second-harmonic  generation 

Second-harmonic  generation  requires  a  strong  incident 
second-subharmonic  field.  As  a  reference  we  consider  an 
incident  power  PSF  of  the  second-subharmonic  field  to 
be  2  W  and  periodical  poling.  We  have  for  the  principal 
squeeze  variances  XSF  «  0.87  and  XPF  «  0.99  for  the 
period  of  poling  Krnl  «  3.552  x  10-3  that  assures  a  perfect 
quasi-phase  matching  (see  Fig.  4b). 

The  nonlinearly  interacting  fields  behave  similarly  as 
in  the  case  of  second-subharmonic  generation.  If  we  in¬ 
troduce  a  linear  corrugation  into  the  second-subharmonic 
field  and  set  the  nonlinear  phase  mismatch  5 pl  1  equal  to 
3.53  x  10-3  (see  Fig.  4b),  values  of  the  principal  squeeze 
variance  XSF  approach  0.6  inside  the  strips  in  the  plane 
spanned  by  variables  Krs  and  5rs  where  the  conditions 
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FIG.  7:  Principal  squeeze  variance  X3F  (a)  and  number  NSp 
of  photons  leaving  the  waveguide  for  the  forward-propagating 
second-subharmonic  field  as  they  depend  on  the  incident 
pump-field  power  Ppf  under  different  conditions:  perfect 
quasi-phase  matching  (<5£;  i  =  0)  without  a  corrugation 
(solid  curve),  quasi-phase  matching  with  Xrnl  =  3.53  x  10-3 
(Sni  i  =  —10.82)  without  a  corrugation  (solid  curve  with  o), 
and  quasi-phase  matching  with  =  3.53  x  10-3  together 
with  a  corrugation  [its  parameters  are  given  by  the  condi¬ 
tions  in  Eqs.  (43)  (m  =  1)  and  (35)]  (solid  curve  with  o); 
PPF  =  2  W,  PSF  =  1  x  10”10  W,  arg(A'p)  =  tt/2;  values  of 
the  other  parameters  are  the  same  as  in  Fig.  4. 


from  Eqs.  (38)  and  (40)  are  fulfilled.  The  number  NSp  of 
second-subharmonic  photons  at  the  output  of  the  waveg¬ 
uide  (together  with  the  output  power)  decreases  by  an 
order  of  magnitude  inside  these  strips  compared  to  the 
incident  power  as  a  consequence  of  transfer  of  energy 
from  the  strong  second-subharmonic  field  into  the  pump 
field  (due  to  an  efficient  nonlinear  interaction)  and  also 
transfer  of  energy  into  the  backward-propagating  field 
(due  to  scattering)  is  considerable.  Despite  this  the 
squeezed  second-subharmonic  field  remains  very  strong, 
it  contains  about  106  photons  inside  the  waveguide.  An 
appropriate  choice  of  the  value  of  nonlinear  phase  mis¬ 
match  Spl  1  is  governed  by  the  same  rules  as  mentioned 
for  second-subharmonic  generation.  The  pump  field  that 
is  only  weakly  squeezed  for  perfect  quasi-phase  matching 
can  reach  values  of  the  principal  squeeze  variance  XPF 
around  0.8  assuming  corrugation  with  parameters  given 
by  Eqs.  (38)  and  (40).  The  pump  field  gets  a  consider¬ 
able  amount  of  energy  from  the  second-subharmonic  field 
and  so  typical  values  of  the  number  NPf  of  pump  pho¬ 
tons  leaving  the  waveguide  can  reach  106;  i.e.  the  output 
power  is  of  the  order  of  10“ 1  W  (energy  of  one  pump  pho¬ 
ton  inside  the  waveguide  of  thickness  t  =  5  x  10“ '  m  and 
length  L  =  lx  10-3  m  corresponds  to  the  output  power 
of  5.13  x  10-8  W).  A  detailed  behavior  of  the  pump  field 
with  respect  to  the  parameters  of  the  linear  corrugation 
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FIG.  8:  Contour  plots  of  principal  squeeze  variance  XPF  (a) 
and  number  NPF  of  photons  leaving  the  waveguide  (b)  for 
the  forward-propagating  pump  field  as  they  depend  on  pa¬ 
rameters  \Kg  \  and  8ra  of  the  linear  corrugation  in  the  second- 
subharmonic  field;  more-less  equidistant  curves  in  graphs  indi¬ 
cate  positions  of  transmission  peaks  [Eq.  (40)],  the  last  curve 
going  up  is  given  by  Eq.  (38);  Arnl  =  3.53  x  10-3  (5^,4  = 
-10.82),  PPF  =  1  x  lO”10  W,  PSF  =  2  W,  arg(AJ)  =  tt/2; 
values  of  the  other  parameters  are  the  same  as  in  Fig.  4.  As  for 
triangles  in  lower  right  corners  of  both  graphs  with  formally 
zero  values,  see  caption  to  Fig.  5. 


can  be  understood  from  graphs  in  Fig.  8. 

A  linear  corrugation  can  also  be  put  into  the  pump  held 
and  we  arrive  at  qualitatively  the  same  results  as  if  the 
corrugation  is  present  in  the  second-subharmonic  held. 
For  the  value  of  nonlinear  phase  mismatch  5^  1  =  —10.82 
used  above,  values  of  the  principal  squeeze  variance  XSp 
can  reach  even  0.3.  On  the  other  hand,  values  of  the 
principal  squeeze  variance  XPF  lie  above  0.9. 

At  the  end  of  discussion  we  make  a  note  about  the 
influence  of  some  phases  with  respect  to  squeezed-light 
generation.  It  can  be  easily  shown  that  any  solution  of 


Eqs.  (8)  depends  only  on  phase  ip  =  arg(Ap)  —  2  arg (Ks). 
In  our  numerical  investigations,  we  did  not  observe  any 
dependence  of  numbers  N  of  photons  as  well  as  principle 
squeeze  variances  A  on  the  phase  ip.  On  the  other  hand, 
these  quantities  depend  weakly  on  phases  arg(ASjr)  and 
arg( ApF)  of  the  incident  fields.  However,  this  depen¬ 
dence  is  very  weak  under  the  conditions  where  a  strongly 
squeezed  light  is  generated. 


V.  MODEL  FOR  A  MULTI-MODE  NONLINEAR 
INTERACTION 


We  generalize  the  model  presented  for  cw  pumping  in 
Sec.  II  to  multi-mode  fields  using  quantum  theory.  Be¬ 
cause  of  the  necessity  to  include  also  inter-mode  disper¬ 
sion,  we  decompose  the  electric-field  operator  amplitudes 
Ea  (a  =  p,  s )  using  mode  operator  amplitudes  da  in  the 
Heisenberg  picture  [5,  8]: 

/•OO 

E  a(x,y,z,t)  =  /  du>a'Ea(x,y,Z,UJa)exp(-ioJat), 

Jo 

(48) 

'Ea(x,y,z,uJa)  =  i[aaF(z,uJa)ea(x,y,u>a ) 

+  aaB(z,ua)ea(x,y,u>a)  -  H.c.] ,  a  =  p,s,  (49) 


where  aaF(z,wa)  [a0j3  (z,  u>a)\  denotes  annihilation  oper¬ 
ator  of  mode  with  frequency  u>a  in  field  a  propagating 
forward  [backward] .  We  note  that  mode  functions  ea  de¬ 
rived  in  Appendix  A  as  well  as  propagation  constants  (3a 
depend  on  the  frequency  ioa  in  the  multi-mode  model. 

Evolution  of  quantum  optical  fields  inside  the  studied 
waveguide  is  described  using  the  following  momentum 
operator  G  [8,  20]: 


/•OO 

G(z)  =  /  dujahpab(oja)aib(z,ua)aab(z,uja) 

a=p,s  b=F,B  v'° 

/•OO 

^  /  du>a^Ka(LOa)alF(z,LOa)daB(z,UJa)  +  H.C. 

a=p,s  ^  0 

nOO  nOQ 

^2 i  dujs  /  du>sKnitq(u)s,us) 

I _ TP  D  J  0  JO 


b—F,B 


X  a\  {z,ws)a\  {z,Js)aPb{z,ua  +J8)  +  H.c.  .  (50) 


The  Heisenberg  equations  written  in  Eq.  (13)  and  corre¬ 
sponding  to  the  momentum  operator  G  in  Eq.  (50)  take 
the  form: 


dciSF{z,u>s) 

dz 


d&sB{z,  u>s) 
dz 


=  if3s{ojs)aSF{z,uis)  +iKs(u>s)aSB(z,us ) 
dujsKnUq(ojs,us)dPF  (z,  us  +  vs)afSF  (z,  ua), 

=  - if3s(u>s)aSB(z,us )  -  iK*(u>s)dSF(z,ojs) 
duj'sKnhq(ujs,u)'s)aPB  (z,  ujs  +  u>'s)dlB  (z,  w'), 
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daPF(z,up) 

dz 


—  (z,  cjp)  H-  iKp{iOp)dpB  (z,  uJp) 


/»oo 

2  /  dusKnqq{uS)Up  Ua)ciSF  (z,  Us)dSF  (.2,  Up  U^s)i 

J  0 


dcips  (  A  Up) 

dz 


—  iPp(^p)®PB  (A  ^p)  -)(Up)dpF(^Z,Up) 


/*oo 

T2  /  dus  K r)jq(u  s .  Up  us)dSB  (z,  Us)dSB  (.2,  UJp  ^s)- 

•lo 

(51) 


Linear  coupling  constants  A's ,  7\p  and  nonlinear  coupling 
constants  Kni,q  ( <7  =  0,  ±1)  are  given  in  Eqs.  (11)  and 
(12)  and  are  frequency  dependent.  We  note  that  the  non¬ 
linear  coupling  constants  Knj  q  have  the  following  sym¬ 
metry:  Knhq{us,us)  =  Kn^q(us,us). 

Creation  and  annihilation  operator  amplitudes  of  the 
incident  fields  are  assumed  to  fulfill  boson  commutation 
relations,  i.e. 

[^fljr(0i^a)i®fl'  (0?^a7)]  —  3a,a'3(oJa  0Ja'), 

[®aB  (1)  ^a)j  (f  i  ^a7)]  —  ^a,a7^(^a  0Ja>), 

a  =  p,s  (52) 


and  the  remaining  commutators  are  zero.  It  can  be  shown 
[49]  for  quadratic  momentum  operators  G  that  also  the 
output  operator  amplitudes  have  to  obey  boson  commu¬ 
tation  relations: 


[haj?  (L,  Ua) ,  d  a,  (L,UJa')]  da,a'5{u>a  0Ja')j 

[ha^  (O5  OJa),  da,  (0,UJa')]  —  da^a'6(u)a  Ua')j 

a  =  p,s  (53) 


and  commutators  not  mentioned  in  Eqs.  (53)  are  zero. 

The  nonlinear  operator  equations  written  in  Eqs.  (51) 
have  one  integral  of  motion  that  expresses  the  conserva¬ 
tion  law  of  the  overall  energy: 


(54) 


When  solving  evolution  of  the  nonlinearly  interact¬ 
ing  fields  as  described  in  Eqs.  (51)  electric-field  op¬ 
erator  amplitudes  ASf ,  ASB,  APf,  and  APB  in  the 
interaction  representation  can  be  conveniently  defined 
(aaF(z,ua)  =  AaF{z,ua)exp[i(3a(ua)z\,  daB(z,ua )  = 
AaB  (z,  ua)  exp[— i[3a{ua)z\,  a  =  p,s)  transforming  the 
equations  in  Eqs.  (51)  into  the  form: 


dASF(z,us) 

dz 


=  iKs(us)  exp[-i5s(ua)z\ASB(z,  us) 


+  4  /  dusKnhq{us,us)exp{i5n^q{ua,ua)z\ 

Jo 

X  ^ Pf  (  A  ^ S  T  ^s)^SF  (A  ^s)  3 
dASB(z,us)  =  _^*(Ws)exp[^a(u;g)2,]_4  sf(z,ws) 


dz 


dAPF  (z,  up) 
dz 


4  /  dusKn^q(ua,us)exp[-iSni,q(ua,us)z} 

Jo 

X  Apb(z,us  +  us)A\b(z,u's), 

=  iKp(up)  exp[-iSp(up)z\ApB  (z,  up) 


roo 

-2  duaK*lq(us,up  -  ua)exp[-i6n^q(ua,up  -  ua)z) 

Jo 


x  ASf  (z,  us)ASf  (z,  Up  -  us) 
dApB  (2,  Up) 
dz 


=  -iKp(up)exp[i6p(up)z\APF(z,Up) 

/»oo 

+  2  /  dusKBlq{us,  Up  -  us)  exp[iSni}q(us,  up  -  us)z] 
Jo 

x  Asb(z,us)Asb(z,up  -  us).  (55) 

Linear  phase  mismatches  Sp,  6S  and  nonlinear  phase  mis¬ 
matches  Sni^q  ( q  =  0,  ±1)  are  defined  as  follows: 

2tt 

iK)  =  2fJa{ua)  — 7—,  a  =  p,s, 

A; 

5ni,q(us,us)  =  (3p{us  +  us)  -  (3s(us)-  Ps{us)  +  q 


27 r 
Anl 

(56) 


It  holds  that  Sn^q(us,us)  =  5ni,q(us,  ua). 


VI.  SECOND-SUBHARMONIC  GENERATION 
WITH  A  LINEAR  CORRUGATION  IN  A 
STRONG  PUMP  FIELD  -  PERTURBATION 
APPROACH 


A  strong  incident  pump  field  being  scattered  on  a  lin¬ 
ear  corrugation  is  considered.  We  also  assume  that  the 
process  of  second-subharmonic  generation  is  weak  and 
the  generated  second-subharmonic  field  remains  at  single 
photon  level.  Because  the  generated  second-subharmonic 
field  is  many  orders  of  magnitude  weaker  than  the  pump 
field,  depletion  as  well  as  quantum  features  of  the  pump 
field  can  be  neglected  and  the  last  two  equations  in 
Eqs.  (51)  can  be  transformed  into  the  equations  for 
strong  amplitudes  aPF  and  aPB : 


daPF  (z,  Up) 
dz 

daPB  (z,  Up) 
dz 


ifJp(up)aPF(z,Up)  +  iKp(up)aPB(z,Up), 

ifJp(uJp)(lpB  (Z-f  t-Zp)  *Ap  (up)cipF  (z,  Up) . 

(57) 


Solution  for  the  forward-propagating  pump-field  ampli¬ 
tude  aPF  can  be  found  in  the  form: 

aPF  (z,  Up)  =  exp[i(3p(up)z\  exp(—i5pz/2) 
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x  [exp(iApz)B+F(u>p)  +  exp(-iApz)BpF(ajp)]  , 

(58) 

where  functions  Bpp  ( B~F  =  B+* )  are  given  by  boundary 
conditions,  i.e.  by  spectrum  of  the  incident  pump  pulse. 
If  the  central  frequency  u)p  is  tuned  at  a  transmission 
peak  [ApL  =  rmr,  m  =  1,2,...,  <5°  =  (5p(w°)]  we  may 
approximately  write: 

aPF(z,u>p)  =  exp[i/3p(up)z\exp(-idpZ/2) 

x  [exp{iA°pz)B°+(ojp)  +  exp{-iA°pz)B°F(up)]  , 

2A°  +  <5° 

K+>p)  =  4Ao  PApA0,“p)\  (59) 

ApF(0,Wp)  means  the  incident  pump-field  amplitude 
spectrum. 

An  operator  perturbation  solution  of  the  first  equation 
in  Eqs.  (51)  with  the  substituted  pump-field  amplitude 
aPF  as  written  in  Eq.  (59)  is  obtained  in  the  form: 


aSF(L,  us) 

U(u)s,u  J 

V(us,us) 


/  du>sU(us,ujs)aSF(0,u>s) 

Jo 

nOO 

+  /  dw'F(ws,a;')aJF(0,w'), 

Jo 


S(ujs  -  cus)  exp [i/3s(uja)L], 
-i4Kni,q(u>s,  ujs)  exp [if3s(us)L\ 


(60) 

(61) 


L 

x  -sine 
2 


sinc(a:)  =  sin(x)/a:.  Two  different  solutions  are  written 
in  Eqs.  (60 — 62):  one  solution  is  appropriate  under  the 
condition  5ni,q{u>s,uis)—5pl 2+A°  «  0  whereas  the  second 
one  is  valid  when  Snl,q(us,  ujs)  -  5p/2  -  A°  «  0.  We  note 
that  one  of  these  two  conditions  has  to  be  fulfilled  in  order 
to  observe  an  efficient  second-subharmonic  generation  as 
has  been  shown  for  cw  case  in  Sec.  III. 

In  our  first-order  perturbation  approach,  the  Green 
function  U  as  given  in  Eq.  (61)  describes  free-held  evo¬ 
lution  whereas  the  Green  function  V  in  Eq.  (62)  charac¬ 
terizes  properties  of  a  generated  photon  pair.  The  Green 
function  V  is  closely  related  to  the  so-called  two-photon 
amplitude  $  introduced  in  the  description  of  spontaneous 
parametric  down-conversion  [50-52] .  Similarly  as  a  two- 
photon  amplitude  4>  the  Green  function  V  can  be  decom¬ 
posed  into  the  Schmidt  decomposition  that  gives  typical 
eigenmodes  of  the  nonlinear  interaction.  This  decomposi¬ 
tion  can  be  reached  even  analytically  under  the  following 
conditions.  We  first  assume  that  the  nonlinear  compling 
constant  Kni,q  does  not  depend  on  frequencies  and  ex¬ 
press  the  propagation  constants  f3a  in  Taylor  series: 


P°a  =  M), 

Pi, a  =  >  *  =  1>2;  a  =  p,s.  (63) 

d0Ja  u,a=u,S 


The  sine  function  occurring  in  Eq.  (62)  can  be  approxi¬ 
mated  by  a  gaussian  function  along  the  relation: 

sine  cxq  +  oq(ws  +  ios  —  Wp)  +  oii (cns  +  los  —  Wp)2 

+a3(u)3  -  us)2  « 


K* 3  /  '  \2  al 

exp - —  (ws-wj  -  — 

3  5 


■  '  o  ,  ao 

+ - 

ai 


(64) 


where  ao,  on,  0:2  >  and  a3  are  constants  and  «2  can  be 
neglected  for  values  of  parameters  typical  for  the  consid¬ 
ered  waveguide.  An  incident  gaussian  pump  pulse  with 
such  phase  modulation  along  the  spectrum  that  phases 
of  different  frequencies  are  equal  in  the  middle  of  the 
waveguide  is  further  assumed  [34],  i.e. 


APF(0,ujp)  =  APf  exp 


^-i[Pp(uP)  ~  Pp]L' j 


and  tp  denotes  pump-pulse  duration.  Constant  APf  in 
Eq.  (65)  gives  strength  of  the  pump  pulse  and  is  deter¬ 
mined  from  the  incident  pump  power  PPF  along  the  for¬ 
mula  in  Eq.  (45).  We  note  that  /0°°  du>p\APF(0,ujp)\2  = 
|ApF|2/27r.  The  Green  function  V  finally  gets  the  gaus¬ 
sian  form: 


Kl,g  =  KnU  wj,wj).  Parameters  A+  and  A_  of  the 
Green  function  V  are  given  as: 

1  P2,SL 

A2_  24  ’ 


J_  =  (Pl,p~Pl,s)2L2  T2 

A2.  20  +  4  ' 


(69) 
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The  gaussian  function  $  given  in  Eq.  (68)  can  be  de¬ 
composed  into  the  Schmidt  decomposition  [34,  54]: 

OO 

wj  =  ^  Hn(t>n{us)4>n{us).  (70) 

n—0 

Eigenvalues  fin  are  given  as  follows: 


-  Va+a-0  ’ 

(71) 

/  A  i  —  A_  \  2 

(72) 

(A-  •  A  ) 

Eigenmode  functions  <j>n  can  be  expressed  in  terms  of  the 
Hermite  polynomials  Hn: 


&*(“>«)  =  \/2"?x!-  ^exp[-^2K-^)2. 


x  Hn(rs[ws  -  ws]), 


1  -  0‘2 


Ts  = 


AO 

°nl,q 


2  2(/3liP -ft,,)' 


(73) 


(74) 


For  the  considered  waveguide,  A+  <C  A_  and  thus  a 
typical  time  constant  rs  defined  in  Eqs.  (74)  can  be  ap¬ 
proximated  as: 


Ts 


(75) 


Taking  into  account  inter-mode  phase  modulation  as  it 
appears  in  Eq.  (67)  eigenmode  functions  4>n  of  the  non¬ 
linear  interaction  are  given  as  folows: 

4>(u>s)  =  exp  [i(3s  (<jja)L/2]<i>(u}s)  (76) 


FIG.  9:  Absolute  value  of  the  linear  coupling  constant  Kp 
(a)  and  relative  linear  phase  mismatch  Sp  —  5p’°  (b)  for  the 
pump  held  as  they  depend  on  the  relative  frequency  ujv/u}p, 
Sp’°  =  L5p(ujp);  t  —  5  x  10-7  m,  ti  =  5  x  10-8  m. 


contained  in  a  multiplicative  factor  (2Aj]  +  5p)/{ 4A°)  in 
Eq.  (66)  that  can  be  expressed  as  follows: 


2A°  +  a°_l  1  l\K°\*  +  (mn/L)2 
4A£  2  2  Y  ( mn/L )2 


(80) 


The  greatest  value  of  enhancement  factor  M  is  reached 
at  the  first  transmission  peak  (m  =1)  and  the  greater 
the  value  of  linear  coupling  constant  I\p  the  greater  the 
value  of  enhancement  factor  M. .  On  the  other  hand  there 
is  no  enhancement  in  the  limit  m  — >  oo,  i.e.  the  effect  of 
pump-field  localization  is  lost. 


and  the  function  <f>  written  in  Eq.  (67)  attains  its  Schmidt 

decomposition  [34]:  VII.  SECOND-SUBHARMONIC  GENERATION 

WITH  A  LINEAR  CORRUGATION  IN  A 
STRONG  PUMP  FIELD  -  NON-PERTURBATION 
APPROACH 

n—0 


OO 

$(ws,  w'  )  =  ^2  Hn<!>n{Us)<t>*n{Us).  (77) 


The  Schmidt  decomposition  of  function  $  written  in 
Eq.  (77)  suggets  to  introduce  new  discrete  sets  of  annihi¬ 
lation  and  creation  operator  amplitudes  that  diagonalize 
the  nonlinear  interaction: 

pOO 

aTF,n  =  /  du>s(f)n(vs)asF(  0,u;a),  (78) 

Jo 

POO 

a°sf,n  =  /  <fas<t>n{us)aSF{L,us).  (79) 

Jo 

Because  functions  <f>n  ((/>*)  form  an  orthonormal  basis, 
operators  and  alsnJn  (a°“(„  and  fulfill  the 

usual  boson  commutation  relations. 

Enhancement  of  the  nonlinear  interaction  caused  by 
scattering  of  the  pump  field  on  a  linear  corrugation  is 


Linear  coupling  constants  Kp ,  Ks  and  nonlinear  cou¬ 
pling  constants  Kni,q  depend  only  weakly  on  frequen¬ 
cies  in  a  relatively  broad  range.  This  is  documented  in 
Figs.  9a,  10a,  and  11a  where  the  dimensionless  constants 
[for  their  definitions,  see  Eqs.  (47)]  K%,  Kp,  and  j ,  re¬ 
spectively,  are  plotted  as  functions  of  frequencies  in  the 
range  of  cca  20  nm  for  typical  values  of  the  waveguide 
parameters. 

If  inter-mode  dispersion  is  omitted,  equations  for  op¬ 
erator  amplitudes  Aab(z,u>)  given  in  Eqs.  (55)  can  be 
decoupled  using  Fourier  transform: 

„  1  f00  , 

Aab(z,ua)=—  dTaAab(z,T)exp(iu>aTa).  (81) 
J  —  OO 


16 


FIG.  10:  Absolute  value  of  the  linear  coupling  constant 
Ka  (a)  and  relative  linear  phase  mismatch  5ra  —  <5J’°  (b)  for 
the  second-subharmonic  field  as  functions  of  the  relative  fre¬ 
quency  ujs/ujp,  Sa’°  =  L5s{u>s);  values  of  parameters  are  the 
same  as  in  Fig.  9. 


The  system  of  equations  in  Eqs.  (55)  attains  the  following 
form  after  Fourier  transform: 


dASF(z,r) 

dz 

4 Knl,q  eMiSnl,qz)ApF  (z,  t)A\f  (z,  t), 
dASB(z,r) 


dz 


-4 Kni,q  exp (-iS°^qz)ApB  (z,  t)A\b  (z,  t), 


dAPF(z,r) 


±pFy~',;  =  iK°exp(-id°z)APB(z,T) 


dz 


2Kniqe>A>(-itinl,qZ)ASF(z,T), 


dAPB(z,r) 
dz 


=  -iKp*  exp(iSpz)APF(z,T) 


+  2Kniq  exp {idnUqz)ASB  (z,  r); 


(82) 


K°a  =  Ka(u°a),  S°a  =  Sa(v°a)  (a  =  p,s),  AZ,g  = 
Kni,q(u°,u°),  and  S^l  q  =  <5„;,g(w°,  w°).  A  suitable  range 
of  frequencies  in  which  this  approximation  is  valid  can  be 
guessed  from  the  curves  in  Figs.  9b,  10b,  and  lib  showing 
the  dependence  of  relative  pump-field  linear  dephasing 
Sp  —  Sp°,  second-subharmonic  linear  dephasing  <5£  —  SF,°, 
and  nonlinear  dephasing  STn[  1  —  cTf  -j ,  respectively,  on  fre¬ 
quencies.  The  requirement  to  neglect  inter-mode  dispe- 
sion  limits  the  length  L  of  waveguide  to  smaller  values. 

If  the  incident  pump  field  is  strong,  its  depletion  due 
to  nonlinear  interaction  can  be  omitted.  Solution  of  the 
last  two  equations  in  Eqs.  (82)  considered  as  equations 
for  classical  amplitudes  APf  and  APB  can  then  be  found 
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FIG.  11:  Contour  plots  of  the  absolute  value  of  nonlinear 
coupling  constant  K^i  ±  (a)  and  relative  nonlinear  phase  mis- 
match  5rnlA-  5^°  (b)  as  they  depend  on  relative  frequen¬ 
cies  ws/wp  and  u>s/uip.  Values  are  shown  only  in  the  area 
\u)s  +  0JS  ~  Wp  <  Awp,  where  A uiv  gives  a  maximum  al¬ 
lowed  width  of  the  pump-field  spectrum,  i.e.  the  black  lower 
left  and  upper  right  triangles  give  no  information.  <5^*^  = 
L5„i, i(u>a,  uja);  values  of  parameters  are  the  same  as  in  Fig.  9. 


in  a  simple  analytical  form  for  transmission  peaks: 

APF{z,Tp)  =  exp(-iS°z/2) 

x  [exp(iA°2)S°+(rp) +  exp(-iA°2)S°;(Tp)]  , 

K+Atp)  =  MAPf(0,tp).  (83) 

Symbol  Apf(0,tp)  denotes  the  Fourier  transform  of  inci¬ 
dent  pump-field  amplitude  spectrum  APf(0,lup)  and  the 
enhancement  factor  A4  is  defined  in  Eq.  (80). 

The  first  operator  equation  in  Eqs.  (82)  can  be  rewrit¬ 
ten  for  two  resonant  conditions  indicated  in  the  expres¬ 
sion  for  amplitude  APf(z,t )  in  Eq.  (83)  as  follows: 

dASF  (z,  t)  o  o0 ±/_\ 

- ~t~ -  =  ^Knl,qBPF(.T) 
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x  exp[i(<S“i)9  -  <S°/2  ±  A°p)z}A\f{z,t).  (84) 

Substitution  ASf(z,t)  =  ASf(z,t)  exp[i(5°j>g/2  - £°/4± 
Ap/2)^],  Atp(z,r)  =  -4iF(z,T)exp[— i(6°lq/2  -  S°/A± 
A°/2)s]  in  Eq.  (84)  leads  to  differential  equations  with 
constant  coefficients  for  operator  amplitudes  ASf  and 
.4jF  that  allow  us  to  find  solution  of  Eq.  (84): 


Asf(L,t)  =  U(t)Asf(0,t)  +  V(t)A\  (0,t),  (85) 


U(t)  =  ^  exp  (  i 


1  - 


+  1 


2A±(r) 


0±L\ 

“2“  J 

expfA*  (r)L] 


2A±(r)  J 


\n  \  ^±(T) 

^(T)  =  nOU  eXP  I  » 


2A±(r) 


)  exp[— A=t(r)I/] 
0±L\ 


(86) 


2  1 


x  [exp[A±(r)i]  —  exp[— A±(r)L]]  .  (87) 

Phase  mismatches  Q  ,  nonlinear  coupling  constants  /C  , 
and  eigenvalues  are  given  by  the  following  formulas: 

^  =  <&,,  -  «5°/2  ±  A°, 

2A°  +  (5° 

^(r)  =  K°nltq^—^APF(0,r), 

^ p 

A±(r)  =  y/|/C±(r)|2  —  fi±2/4.  (88) 


Inverse  Fourier  transform  of  the  expression  for  oper¬ 
ator  amplitude  ASf  ( L ,  r)  in  Eq.  (85)  provides  the  op¬ 
erator  amplitude  ASf(L,ojs)  and,  after  returning  from 
the  interaction  representation,  the  following  expression 
for  the  output  second-subharmonic  operator  amplitude 
aSF  ( L ,  los)  is  reached: 


Operator  amplitudes  aSFy  defined  in  Eq.  (90)  obey  the 
usual  boson  commutation  relations  instead  of  those  writ¬ 
ten  in  Eq.  (53).  Transformation  written  in  Eq.  (60) 
between  the  input  second-subharmonic  operator  ampli¬ 
tudes  a*™  i,  a*p  i  and  output  second-subharmonic  oper¬ 
ator  amplitudes  a°“b,  a°p*f  can,  after  discretization,  be 
written  in  the  matrix  form: 

a^  =  Ua“+Va“t.  (92) 

Vectors  a*”  and  a°“*  are  composed  of  operator  ampli¬ 
tudes  i  and  a°Fb,  respectively. 

Because  the  input  operator  amplitudes  a*Fji,  a*F)i  as 

well  as  the  output  operator  amplitudes  a°Fti,  a°“*J  ful¬ 
fill  boson  commutation  relations  the  transformation  writ¬ 
ten  in  Eq.  (92)  represents  the  Bogoljubov  transformation. 
This  means  that  matrices  U  and  V  can  be  decomposed 
using  the  Bloch — Messiah  reduction  [33,  34]: 

U  =  XAtfY*, 

V  =  XAvYt,  (93) 

where  f  means  hermitian  conjugated  matrix  and  T  stands 
for  matrix  transposition.  Matrices  A jj  and  Ay  are  diag¬ 
onal  and  contain  nonnegative  eigenvalues  of  the  decom¬ 
position.  Matrix  Y  (X)  in  Eq.  (93)  contains  the  right 
(left)  eigenvectors. 

Eigenvectors  defined  by  the  Bloch-Messiah  reduction 
give  typical  modes  of  the  nonlinear  interaction  and  rep¬ 
resent  a  discrete  form  of  eigenmode  functions  (j)n  found 
in  Eq.  (76)  in  analytical  approach  of  Sec.  VI. 

Coefficents  Bn  and  Cn  of  the  generalized  superposition 
of  signal  and  noise  [8]  valid  for  the  nth  eigenmode  and  de¬ 
fined  in  Eqs.  (24)  can  be  expressed  in  terms  of  eigenvalues 
of  the  decomposition  (an  incident  second-subharmonic 
field  in  vacuum  state  is  assumed): 


clSf(L,los)  =  exp  [i(3(w8)L\ 


du>sU(us  -  u!s)  ASf  (0,  wj 


IV  o 


/»oo 

+  /  dusV(us  +  ws)A^f(0,o;s) 
Jo 


(89) 


Functions  U (o>)  and  V (u>)  are  given  by  inverse  Fourier 
transform  of  the  expressions  occurring  in  Eqs.  (86)  and 
(87)  and  can  be  found  numerically. 

The  general  expression  for  continuous  output  operator 
amplitudes  aSF{L,u>s )  as  it  occurs  in  Eq.  (60)  has  to  be 
discretized  in  order  to  investigate  the  model  numerically. 
We  introduce  discrete  mode  operator  amplitudes  aSFti 
that  are  distant  by  A u>,  i.e. 


dSF,i(z)  =  V Alo  aSF(z,uj°  +  iAoj).  (90) 


The  Green  functions  U  and  V  are  then  transformed  into 
matrices  U  and  V : 


U ij  =  Aw  (7(w°  +  iAw,  +  jAw), 

=  Aw  F(wJ  +  i Aw,  uj  +  jAu) .  (91) 


Bn  =  (Av)L  > 

Cn  =  (Au)nn{ Ay)nn.  (94) 

The  second-subharmonic  principal  squeeze  variance 
A SFin  of  eigenmode  n  is  then  determined  along  the  ex¬ 
pression: 

XSF,n  =  1  +  2  (Av)nn  [(Av)nn  -  (A u)nn)  ■  (95) 

The  number  ASFj„  of  photons  leaving  the  waveg¬ 
uide  in  eigenmode  n  is  given  by  a  simple  formula 
(Nsf  =  {A\fASf),  averaging  is  over  an  incident  second- 
subharmonic  vacuum  state): 

Nsr,n  =  (AV)ln  ■  (96) 

Importance  of  the  Bloch-Messiah  reduction  in  investi¬ 
gations  of  squeezing  is  emphasized  by  the  fact  that  an 
eigenmode  with  the  lowest  value  of  the  principal  squeeze 
variance  A SFy  represents  the  solution  of  the  optimiza¬ 
tion  problem  for  a  suitable  spectral  profile  of  a  mode 
that  gives  the  best  possible  amount  of  squeezing  (see  Ap¬ 
pendix  B). 
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FIG.  12:  Enhancement  factor  A4  as  a  function  of  the  depth  ti 
of  linear  corrugation  is  plotted  for  the  first  live  transmission 
peaks:  m  =  1  (solid  curve),  2  (solid  curve  with  o),  3  (solid 
curve  with  *),  4  (solid  curve  with  A),  and  5  (solid  curve  with 
o);  t  =  5  x  10' '  nr. 


VIII.  SQUEEZED-LIGHT  GENERATION  - 
NUMERICAL  ANALYSIS  FOR  PULSED 
PUMPING 


FIG.  13:  Contour  plot  of  the  principal  squeeze  variance  XSF,n 
as  a  function  of  incident  pump  power  PVF  and  mode  number 
n;  tp  =  1  x  10~13  s,  values  of  the  other  parameters  are  the 
same  as  in  Fig.  12. 


Benefit  of  a  linear  corrugation  present  in  the  pump  field 
in  squeezed-light  generation  can  be  roughly  judged  using 
the  enhancement  factor  M  defined  in  Eq.  (80).  Values 
of  the  enhancement  factor  M  around  4  can  be  reached 
in  the  studied  waveguide  for  the  first  transmission  peak, 
as  documented  in  Fig.  12. 

A  detailed  analysis  of  squeezed-light  generation  in 
pulsed  regime  is  given  for  fixed  values  of  the  depths 
of  waveguide  (t  =  5  x  10“ 7  m)  and  linear  corrugation 
(ti  —  5  x  10~8  m).  The  value  of  enhancement  factor  A4 
equals  cca  1.5  for  the  first  transmission  peak  (see  Fig.  12). 
A  suitable  value  of  the  pump-field  linear  phase  mismatch 
is  given  along  the  formula  in  Eq.  (41)  knowing  the 
value  of  pump- field  linear  coupling  constant  Kp.  It  de¬ 
pends  on  the  number  m  counting  transmission  peaks. 
The  condition  in  Eq.  (35)  then  provides  an  optimum 
value  for  the  nonlinear  phase  mismatch  5^i  q: 

jo  l(S0)2 

=  i  ±\/Mr- W-  (97) 

In  Eqs.  (41)  and  (97)  there  exist  four  possible  combi¬ 
nations  of  signs.  As  for  squeezing,  condition  with  both 
signs  +  gives  the  same  result  as  that  with  both  signs  -. 
Similarly,  conditions  with  one  sign  +  and  one  sign  -  lead 
to  the  same  amount  of  squeezing.  In  our  case,  the  con¬ 
ditions  with  one  sign  +  and  the  other  sign  -  give  better 
values  of  the  principal  squeeze  variance  XSF  and  that  is 
why  we  consider  them  in  the  following  discussion. 

Gaussian  spectrum  Apf(0,lop)  of  the  incident  pump 
pulse  is  assumed.  Its  Fourier  transform  takes  the  form: 


A>F(°>'r)  = 


1  V2 

V^Tp 


exp 


t21 

"(1 1  i(ip )  2 
Tp  ■ 


(98) 


where  tp  means  pump-pulse  duration  and  ap  stands  for 
chirp  parameter. 


In  order  to  observe  squeezing  in  pulsed  regime  the  right 
choice  of  eigenmodes  is  crutial.  It  occurs  that  there  typ¬ 
ically  exist  several  or  maximally  several  tens  of  eigen¬ 
modes  that  can  be  considerably  squeezed.  Values  of  the 
principal  squeeze  variances  XSF,n  of  these  modes  decrease 
with  the  increasing  incident  pump- field  power  PPF,  as 
shown  in  Fig.  13  for  100-fs  pump  pulse.  Measurement 
of  squeezing  in  the  lowest  modes  can  be  reached  using 
a  suitable  profile  of  a  local  oscillator  in  homodyne  de¬ 
tection.  Eigenmode  functions  (j>n  of  the  first  four  modes 
are  depicted  in  Fig.  14.  Amplitudes  of  the  functions  <f>n 
with  even  numbers  n  have  maximum  at  the  central  fre¬ 
quency  a;!?  whereas  those  with  odd  numbers  n  are  zero. 
The  larger  the  mode  number  n  the  faster  the  oscillations 
along  the  frequency  axis  u)s.  The  eigenmode  functions  <j>n 
plotted  in  Fig.  14  are  not  damped  as  a  consequence  of  ne¬ 
glecting  inter-mode  dispersion.  As  the  analytical  results 
given  in  Eqs.  (73)  and  (76)  suggest  the  inclusion  of  inter¬ 
mode  dispersion  leads  to  damped  eigenmode  functions  <f>n 
with  a  typical  time  constant  rs  written  in  Eq.  (74).  We 
note  that  the  Bloch-Messiah  reduction  given  in  Eq.  (93) 
can  be  done  such  that  the  left  and  right  eigenvectors 
(contained  in  matrices  X  and  Y)  coincide. 

Shortening  of  the  pump  pulse  keeping  its  energy  fixed 
results  in  lower  values  of  the  principal  squeeze  variances 
XSF,n  (see  Fig.  15).  This  reflects  the  fact  that  shortening 
of  the  pump  pulse  leads  to  concentration  of  its  energy 
into  a  smaller  area  inside  the  nonlinear  medium  and  this 
increases  efficiency  of  the  nonlinear  process.  Lower  val¬ 
ues  of  the  principal  squeeze  variances  XSF}n  then  natu¬ 
rally  occur.  Alternatively  we  may  argue  that  the  shorter 
the  pump-pulse  duration  tp  the  broader  the  pump-pulse 
spectrum  APF{ 0,u>p)  and  also  the  better  the  ability  to 
excite  the  lowest  eigenvectors  </>„  having  the  lowest  val¬ 
ues  of  XSFtn.  In  cw  limit  (rp  — ►  oo),  the  value  of  ASF,o 
belonging  to  the  lowest  eigenmode  equals  cca  0.9  when 
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FIG.  14:  Amplitude  squared  \4>n\2  (a)  and  phase  arg ((/>n)  (b) 
of  the  first  four  eigenmode  functions  </>n  as  they  depend  on 
the  normalized  frequency  t os/cjp;  n=l  (solid  curve  with  o),  2 
(solid  curve  with  *),  3  (solid  curve  with  A),  and  4  (solid  curve 
with  o);  rp  =  1  x  10-13  s,  PPF  =  2  W,  values  of  the  other 
parameters  are  the  same  as  in  Fig.  12. 
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FIG.  16:  Contour  plots  of  the  principal  squeeze  variance  ASF)o 
(a)  and  number  NSf,o  of  photons  leaving  the  waveguide  (b) 
as  they  depend  on  incident  pump  puwer  PpF  and  pump-pulse 
duration  tp  for  n  =  0.  Logarithmic  scale  is  used  on  the  a  axis 
in  (b).  Values  of  parameters  are  the  same  as  in  Fig.  12. 


FIG.  15:  Contour  plot  of  the  principal  squeeze  variance  A SF,n 
as  it  depends  on  the  pump-pulse  duration  tp  and  mode  num¬ 
ber  n;  PPF  =  2  W,  values  of  the  other  parameters  are  the 
same  as  in  Fig.  12. 

pumped  by  the  power  of  2  W.  On  the  other  hand,  cw 
model  discussed  in  Sec.  IV  gives  the  value  of  XSF  cca  0.3, 
i.e.  cw  model  predicts  higher  squeezing  of  the  second- 
subharmonic  field.  The  main  difference  between  cw  and 
pulsed  models  is  that  energy  converted  from  the  pump 
field  goes  into  many  modes  of  the  second-subharmonic 
field  in  the  pulsed  model  whereas  there  is  only  one  (effec¬ 
tive)  mode  in  the  second-subharmonic  field  in  cw  model. 
Thus  many  modes  in  the  pulsed  model  compete  and  as 
a  consequence  they  cannot  reach  such  low  values  of  the 


principal  squeeze  variance  XSF  as  those  provided  by  cw 
model. 

Only  the  lowest  eigenmode  can  be  measured  in  op¬ 
timum  homodyne  detection  in  which  a  local  oscillator 
has  shape  of  the  eigenmode  function  fo.  Attainable  val¬ 
ues  of  the  principal  squeeze  variance  ASf,o  for  a  broad 
range  of  values  of  pump-pulse  duration  rp  and  incident 
pump-pulse  power  PPF  are  shown  in  Fig.  16a.  The  num¬ 
ber  NSf_q  of  photons  leaving  the  waveguide  depicted 
in  Fig.  16b  clearly  shows,  after  comparison  with  the 
graph  in  Fig.  16a,  that  the  better  the  conversion  of  en¬ 
ergy  into  the  lowest  eigenmode  the  better  the  squeez¬ 
ing  in  this  mode.  For  example,  an  incident  pump  pulse 
having  100  fs  and  power  of  2  W  generates  a  second- 
subharmonic  field  with  the  lowest  eigenmode  reaching 
the  value  ASFiq  =  0.65. 
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The  obtained  results  show  that  the  waveguide  made  of 
LiNbC>3  can  produce  pulsed  squeezed  light  being  emit¬ 
ted  in  the  process  of  second-subharmonic  generation.  A 
linear  periodic  corrugation  with  appropriate  parameters 
and  sensitive  to  a  pump  held  enhances  considerably  the 
nonlinear  process. 

The  analyzed  waveguide  can  serve  also  as  a  source  of 
light  with  sub-Poissonian  photon-number  statistics.  Ori¬ 
gin  of  these  nonclassical  states  having  suppressed  fluc¬ 
tuations  in  their  photon  numbers  lies  in  nonlinear  inter¬ 
action,  similarly  as  in  the  case  of  squeezed-light  genera¬ 
tion.  Suitable  conditions  for  generation  of  light  with  sub- 
Poissonian  photon-number  statistics  are  similar  to  those 
needed  for  observation  of  squeezed  light  and  analyzed  in 
detail  in  the  report.  This  similarity  has  been  elucidated 
in  detail  in  [45,  46]  for  three-mode  interaction.  It  allows 
to  use  directly  the  results  obtained  for  squeezed  light  also 
in  the  analysis  of  suitable  conditions  for  the  generation 
of  light  with  sub-Poissonian  photon-number  statistics. 


IX.  CONCLUSIONS 


We  have  shown  that  an  additional  scattering  of  two 
nonlinearly  interacting  optical  fields  caused  by  a  small 
linear  periodic  corrugation  at  the  surface  of  a  waveguide 
can  lead  to  an  enhancement  of  the  nonlinear  process  thus 
resulting  in  higher  generation  rates  and  better  values  of 
squeezing.  To  observe  this  effect  the  natural  nonlinear 
phase  mismatch  has  to  be  nonzero  in  order  to  match  with 
periodic  oscillations  caused  by  scattering  at  the  corruga¬ 
tion.  Under  these  conditions,  better  values  of  squeezing 
as  well  as  higher  values  of  intensities  of  both  interacting 
fields  (pump  and  second-subharmonic)  can  be  reached  in 
comparison  with  a  perfectly  quasi-phase-matched  waveg¬ 
uide.  A  linear  corrugation  can  be  designed  to  match  ei¬ 
ther  the  pump  or  the  second-subharmonic  field,  or  even 
both  of  them.  Origin  of  the  enhancement  of  nonlinear 
interaction  lies  in  constructive  interference  of  scattered 
light  leading  to  spatial  localization.  Optimum  conditions 
for  this  enhancement  have  been  found  approximately  us¬ 
ing  analytical  approach  and  confirmed  numerically.  In 
practise,  the  best  results  are  reached  with  the  deepest 
possible  corrugation.  Periods  of  linear  corrugation  and 
periodical  poling  are  then  determined  from  the  model. 
In  pulsed  regime  a  suitable  choice  of  eigenmodes  in  the 
second-subharmonic  field  is  crutial  to  reach  high  levels 
of  squeezing.  The  shorter  the  pump  pulse  the  better  the 
values  of  squeezing. 

The  obtained  results  have  shown  that  nonlinear  planar 
waveguides  with  periodically  corrugated  surfaces  repre¬ 
sent  a  promising  source  of  squeezed  light  for  integrated 
optoelectronics  under  the  conditions  excluding  nonlinear 
medium  embedded  inside  a  cavity. 


APPENDIX  A:  MODES  OF  AN  ANISOTROPIC 
WAVEGUIDE 


Mode  of  the  considered  waveguide  [53]  depicted  in 
Fig.  1  is  given  as  a  solution  of  the  wave  equation  written 
in  Eq.  (7).  The  waveguide  is  made  of  LiNbC>3  crystal  us¬ 
ing  the  method  of  proton  exchange.  The  crystallographic 
z  axis  coincides  with  the  x  axis  of  the  coordinate  system 
(see  Fig.  1).  Ordinary  (nS)0)  and  extraordinary  (ns>e)  in¬ 
dices  of  refraction  of  LiNbC>3  valid  for  the  substrate  and 
used  in  calculations  are  given  as: 

n2s,a  =  Aa+A2_aC<  A2,  d  =  0,e, 

Aa  =  4.91300,  Ba  =  0.118717, 

C0  =  0.045932,  D0  =  0.0278, 

Ae  =  4.57906,  Be  =  0.099318, 

Ce  =  0.042286,  De  =  0.0224;  (Al) 

wavelength  A  is  in  /tm.  After  proton  exchange,  ordinary 
(nw>0)  and  extraordinary  (nWif,)  indices  of  refraction  of 
LiNbOa  characterizing  the  waveguide  are  reached: 

1T'W,o  —  ^s,o 
Tlw,e  —  ^s,e  A 

(Sn)~  =  Ai  +  ^  —  DiX2-, 

A-,  =  0.007596,  Bi  =  0.001129, 

Ci  =  0.116926,  D1  =  -0.0003126.  (A2) 

We  assume  that  air  is  present  above  the  waveguide,  i.e.: 

nu  =  1.  (A3) 


Because  only  the  extraordinary  index  of  refraction 
of  LiNbC>3  increases  during  proton  exchange,  only  TM 
waves  can  be  guided.  For  this  reason,  instead  of  solving 
Eq.  (7)  for  x  and  z  components  of  the  electric- field  mode 
functions  ea,  we  solve  the  following  equation  for  the  only 
nonzero  y  component  of  the  magnetic-field  mode  func¬ 
tions  ha  (x)  (fields  are  assumed  to  be  homogeneous  along 
the  y  axis)  [41]: 


d2  ha(x) 
dx2 

(A4) 


ezz(x,L>a)  „2  ,  Cz(x,0Ja)0J2 


p  (x,  UJa) 


P 


ha(x)  =  0. 


Solution  of  Eq.  (A4)  for  the  y  component  of  magnetic- 
field  mode  function  ha(x)  (a  =  p,s)  can  be  written  as: 


[hQ(x)]y  =  —Ca-A  exp(— gx),  x  >  0; 
Qa 


=  C„ 


- -  cos (hax)  +  sin (hax) 

Qa 


cos (hat)  +  sin (hat) 
Qa 


=  -c, 

x  exp(pax),  x  <  — t 


,  0  <  x  <  — t; 
exp  (p  at) 

(A5) 
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where  Ca  denotes  a  normalization  constant.  We  have  the 
following  expressions  for  coefficients  ha,  qa,  pai  pa,  and 
qa  for  the  considered  orientation  of  LiNbCU: 


Solution  written  in  Eq.  (A5)  holds  provided  that  the  fol¬ 
lowing  dispersion  relation  giving  the  propagation  con¬ 
stant  (3a  as  a  function  of  the  frequency  ua  is  valid: 

tan  (hat)  =  — - (A7) 

K  -  paqa 

We  note  that  possible  solutions  of  Eq.  (A7)  for  the  prop¬ 
agation  constant  (3a  lie  in  the  interval  ?rSie(wa)wa/c  < 
Pa  <  nWte(ua)u;a/c. 

Components  of  the  electric-field  mode  function  ea(x) 
can  be  derived  from  the  magnetic-field  mode  functions 
ha(x)  along  the  relations: 


[ea(a:)]x 

[ea(2;)]y 

\ea{x)]z 


Pa 


^a) 

0, 


ixae0ezz{x1uja)  dx 


[ha(x)]y) 


dfh^x)],,  (A8) 


The  normalization  constants  Ca  occurring  in  Eqs.  (A5) 
are  determined  from  the  condition  that  the  mode  func¬ 
tions  ea(x)  describe  one  photon  with  energy  fuxa  inside 
the  waveguide  (of  length  L  and  thickness  Ay): 


2e0A yL  dx  exx(x,ua)  |[ea(x)]3 


z(x,ua)  |[ea(x)]2 


=  Uuj„. 


(A9) 


Corrugation  on  the  surface  causes  periodic  changes  of 
values  of  permittivity  e  for  x  £  (0,  —  ti)  (see  Fig.  1)  and 
we  have  e±  1  =  i(e— 1) /(ne)  in  this  case  using  Eqs.  (2)  and 
(3).  If  the  waveguide  is  periodically  poled  d±i  =  —  2i/-7rd 
in  Eq.  (5)  and  the  remaining  coefficients  may  be  omitted. 
Using  the  electric-field  mode  functions  ea  determined  in 
Eqs.  (A8),  linear  ( Ks ,  Kp)  and  nonlinear  ( Kni  0 ,  Kni  ±i) 
coupling  constants  defined  in  Eqs.  (11)  and  (12)  can  be 
rearranged  into  the  form: 


Kn 


2nc2Pa 


nl,e(Ua)  ~  1 
<i,ek) 


dx  |[ea(x)]x|2 


d^nl.O 


Knl,±  1 


Nonzero  coefficients  of  the  nonlinear  tensor  d  of 
LiNbCU  used  in  calculations  are  the  following: 


=  dZyy  —  dyyZ  —  3.1  X  10  ^mV  \ 

=  dxzz  —  dzzx  —  dyyX  —  5.87  x  10  mV  , 

=  41.05  x  10-12mV_1.  (A13) 


APPENDIX  B:  AN  OPTIMUM  MODE  FOR 
PULSED  SQUEEZED  LIGHT 

We  look  for  a  suitable  linear  combination  of  output 
operator  amplitudes  that  minimizes  the  value  of 

principal  squeeze  variance  XSF  [54]: 

=  (B1) 

i 

Coefficients  t,  fulfill  the  normalization  condition 
\ti\2  =  1.  Using  Eqs.  (24)  and  (26)  the  principal 
squeeze  variance  XSF  as  a  function  of  U  and  t*  can  be 
expressed  in  the  following  form: 

A8F(t,t*)  =  l  +  2ttW1t-2|tW2t|.  (B2) 

Vector  t  occurring  in  Eq.  (B2)  is  composed  of  coefficients 
ti .  Matrices  Wj  and  W2  are  defined  as: 

Wi  =  V*VT, 

W2  =  UVT,  (B3) 

where  the  matrices  U  and  V  are  given  in  Eq.  (91). 

Instead  of  using  coefficients  ti ,  we  define  new  coeffi¬ 
cients  Si  along  the  relation: 

tTWt 

s  =  texp(-iv?/2),  exp(?V)  =  ^Twt^.  (B4) 

Minimum  value  of  the  principal  squeeze  variance  XSF  is 
reached  for  the  values  of  coefficients  Sj  and  s*  that  min¬ 
imize  the  function  /: 

/( s,  s*)  =  1  +  2s^ W is  -  2sW2s  -  Atst  •  s,  (B5) 
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/.t  denotes  a  Lagrange  multiplier  for  the  condition  •  s  = 

1. 

Derivation  of  the  function  /  given  in  Eq.  (B5)  with 
respect  to  coefficients  s*  results  in  the  condition: 

2WiS  =  n  s.  (B6) 

This  condition  can  be  rearranged  using  Eqs.  (B3): 

VV^*  =  |s*.  (B7) 

On  the  other  hand,  the  Bloch-Messiah  reduction  of  ma¬ 
trix  V  in  Eq.  (93)  provides  the  formula 

VV'i’Xj  =  (AvjjjXj  (B8) 

valid  for  itli  eigenvalue  (Av)ii  and  eigenvector  X;.  Com¬ 
parison  of  Eqs.  (B7)  and  (B8)  leads  to  the  conclusion  that 
the  complex  conjugated  vector  s*  coincides  with  the  zth 
eigenvector  X;  and  the  Lagrange  multiplier  /i  fulfills  the 
relation  /z  =  2(Av)fi. 

Because  we  look  for  a  minimum  value  of  the  principal 
squeeze  variance  XSF  the  appropriate  mode  is  character¬ 
ized  by  the  largest  eigenvalue  (Av)i*  [compare  Eq.  (95)]. 


This  mode  has  also  to  fulfill  the  condition  obtained  by 
derivating  the  function  /  written  in  Eq.  (B5)  with  respect 
to  coefficients  s^: 

2sfWi  -  2stW2  -  2W2s  -  /zsf  =  0.  (B9) 

The  use  of  expresions  in  Eqs.  (B3)  and  (B6)  transforms 
Eq.  (B9)  into  the  final  form: 

VUTs  +  sTUVT  =  0.  (BIO) 
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